ASSESSMENT OF A WAKE VORTEX 


FLIGHT TEST PROGRAM 


By Selden B. Spangler, Marnix F. E. Dillenius, 
Richard G. Schwind and Jack N. Nielsen 
Nielsen Engineering & Research, Inc. 


NEAR TR 48 
July 1974 


Prepared under Contract No. NAS2-6973 by 
NIELSEN ENGINEERING & RESEARCH, INC. 
Mountain View, California 


! ‘ ' ' ' ; 'f o± Ames Research Center 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TABLE OF CONTENTS 

Page 

No. 

SUMMARY 1 

INTRODUCTION 1 

LIST OF SYMBOLS 2 

WIND TUNNEL EXPERIENCE ON WAKE VORTEX MEASUREMENTS 5 

NOMINAL FLIGHT VORTEX CHARACTERISTICS 7 

ATMOSPHERIC TURBULENCE MEASUREMENTS 11 

DATA REQUIREMENTS FOR THE FLIGHT TESTS 17 

Airplane Data 18 

Trailing Vortex Data 18 

Meteorological Data 24 

INTERACTION BETWEEN THE PROBE AIRCRAFT AND THE VORTEX WAKE 25 

Steady Flow Disturbance Due to the Aircraft 25 

Aircraft Motion Induced by the Vortex System 26 

Effect of the Aircraft on Vortex Deformation and Stability 28 

FLIGHT TEST PROBE PROCEDURE 30 

EFFECT OF CONDENSATION ON VORTEX CHARACTERISTICS 31 

VORTEX VELOCITY DATA REDUCTION 36 

Methods for Reconstruction of the Vortex Characteristics 37 

Trailing Vortex Crossflow Program 42 

Calculation procedure 43 

Program operation 44 

Program limitations and precautions 44 

Description of input 45 

Description of output 47 

Program listing 49 

Sample Case 57 

CONCLUDING REMARKS 57 

REFERENCES 59 

TABLE I THROUGH IV 61 

FIGURES 1 THROUGH 16 . 64 


iii 


ASSESSMENT OF A WAKE VORTEX FLIGHT TEST PROGRAM 

By Selden B . Spangler, Marnix F. E. Dillenius, 

Richard G. Schwind and Jack N. Nielsen 
Nielsen Engineering & Research, Inc. 

SUMMARY 

An investigation was made of a proposed flight test program by the 
Ames Research Center, NASA to measure the characteristics of wake vortices 
behind a T-33 aircraft. A number of facets of the flight tests were 
examined to define the parameters to be measured, the anticipated vortex 
characteristics, the mutual interference between the probe aircraft and 
the wake, the response of certain instruments to be used in obtaining 
measurements, the effect of condensation on the wake vortices, and methods 
of data reduction. 

On the basis of the investigation, it is recommended that the probe 
aircraft (a Lear jet) make horizontal passes through the vortex wake, 
entering the wake at right angles to its length. Since the generating 
and probe aircraft are. similar in speed, size, and weight, the motions of 
the probe aircraft approaching and passing through the wake are not severe 
and the corrections to the hot-wire data due to aircraft motion should be 
small. In addition to the basic hot-wire measurements of the velocity 
field in the vortices, it is recommended that supporting data be obtained 
to define the environment in which the vortices persist and dissipate. 
Specific recommendations on the type of data and means of measurement are 
given. A data reduction method was developed to reconstruct the velocity 
field and essential features of the vortex, based on knowledge of the 
velocity components along some flight path through the vortex. Finally, 
it is noted that very little information is available on the dissipation 
of vortices, and measurements in this region would be most valuable as an 
aid in understanding the dissipation mechanisms. 

INTRODUCTION 

This report presents the results of a study to assist the Ames 
Research Center staff in defining a flight test program to critically 
evaluate theories for trailing vortex behavior and decay. The objective 



of the study is to examine the various aspects of the proposed flight 
test program to establish requirements for data and to make recommendations 
concerning test procedures, instrumentation, and data reduction. The back- 
ground for this work consists of theoretical analysis on trailing vortices 
done under contract to the Air Force (Contract F44620-70-C-0052) and 
experimental work on measurement of trailing vortex characteristics with 
hot-wire probes in the Ames Research Center 40- by 80-Foot Wind Tunnel 
(Contract NAS2-6719) . 

A number of problems connected with the flight tests were investigated. 
The various theories and analyses of trailing vortices were examined to 
define the parameters that should be measured. Investigations were made 
of atmospheric turbulence, the mutual interference between the probe air- 
craft and the vortex system, effects of condensation on the vortex, the 
response of certain instruments, and methods of data reduction. The 
results of these investigations are described in the various sections of 
this report. In addition, a section is included describing the NEAR 
experience with the 40- by 80-Foot Wind Tunnel program to the extent that 
it could be instructive in planning and executing the flight test program. 

LIST OF SYMBOLS 


b 


c 

c 


P 


c 


wing span 

speed of sound or chord 

specific heat at constant pressure 

wing mean aerodynamic chord 

wing drag coefficient, based on S^ 


C„ rolling moment coefficient, defined as the rolling moment 

divided by qS c 
w 

d distance of a crossflow plane of interest downstream of the 

wing 

e voltage output of a hot-wire probe 

E(k) one-dimensional power spectral density for atmospheric 

turbulent velocity fluctuations along the aircraft 
direction of flight 
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f 

G(k) 


I 

o 

k 

K 

m 


M 

Nu 



Re 


t 

T 


S 


w 


u,v,w 


frequency, Hz 

one-dimensional power spectral density for atmospheric 
turbulent velocity fluctuations normal to the aircraft 
direction of flight 

total enthalpy, defined by equation (11) 

inverse of the wave length, or thermal conductivity 

vortex strength parameter , vr 

mass 

slope of the line normal to the resultant velocity vector V . 
m^ = tan 0^ (figure 8) r 

molecular weight 
Nusselt number 
Prandtl number 

static pressure 

free-stream dynamic pressure, 1/2 pU 2 
radius 

outer radius of vortex eye region, figure 1 

outer radius of vortex logarithmic region, figure 1. 

radius at which K has reached 99 percent of its asymptotic 
value, figure 1 

radius for peak tangential velocity, figure 1 

Reynolds number, based on hot-wire diameter unless otherwise 
indicated 

time 

static temperature 
wing area 

vortex velocities in a cylindrical coordinate system, u radial, 
v tangential, and w axial 
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u'/V'fW'. fluctuating velocity components of atmospheric turbulence, 

u' in direction of aircraft motion, v' lateral to direction 
of aircraft motion and horizontal, and w' vertical 


u ,v ,w steady perturbation velocities induced around the Learjet 
" “ 3 aircraft, u axial, positive forward , v lateral, positive 

cl ci 

to starboard, w^ vertical, positive down 


U 


aircraft flight speed or wind-tunnel speed 


V 


V 


r 


V 

s 


V 


t 


component of the vortex induced perturbation velocity normal 
to the probe flight path for a horizontal probe traverse in 
the crossflow plane 

vector sum of V and V 
n s 

component of the vortex induced perturbation velocity coaxial 
with the probe flight path for a horizontal probe traverse 
in the crossflow plane 

vortex tangential velocity (equal to V ) 


W 


x,y,z 


Z,Y 


a 

a ' 


a 


c 



0 

s 


weight 

axis system with origin at the fuselage nose, x positive 
forward along fuselage centerline, y spanwise positive to 
left, and z positive down 

coordinate system in the crossflow plane of the trailing vortex 
wake, Y positive to the right, Z positive down when looking 
upstream (figure 8) 

aircraft angle of attack 

perturbation angle due to turbulence, equation (6) 

angle between V, and V . Tan a = V /V (figure 8) 
u n c s n 

vortex circulation strength at r Q 

vortex circulation strength at r z 

mean square displacement of the vortex from its mean position 
due to turbulence 

angle between the Z axis and the radius vector to the 
apparent vortex center, measured from the positive Z axis 
in the counterclockwise direction (figure 8) 

angle between the Z axis and the local tangent to the flight 
path, measured from the positive Z axis in the counterclock- 
wise direction. Its tangent is the slope of probe flight path, 
(figure 8) . 
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A 

s 

V 

v 

P 


a 

0) 


a 

ax 

c 

m 

sat 

w 


scale length of turbulence 
viscosity 

kinematic viscosity 
mass density 

rms turbulence velocity, normalized by U 
absolute humidity 

Subscripts 

air 

axial component 
vortex centerline 
mean 

saturation value 
wire or water vapor 

condition at a large radial distance from the vortex 


WIND TUNNEL EXPERIENCE ON WAKE VORTEX MEASUREMENTS 

In January and June of 1972, tests were conducted in the Ames Research 
Center 40- by 80-Foot Wind Tunnel to measure velocity profiles in the 
trailing vortex system behind two wings. The instrumentation used con- 
sisted of a three-wire DISA hot-wire probe . NEAR personnel were intimately 
involved in the design, execution, and data reduction in these tests. The 
experience gained in these tests, insofar as it may have a bearing on the 
proposed flight test program, is briefly described in this section. 

Two rectangular wings were used: an 8-foot span with 18-inch chord 

and a 32-inch span with 6-inch chord (one-third scale) . The wings were 
mounted in the center of the tunnel cross section and angles of attack of 
8° and 12° were tested. Both wings contained smoke generators con- 
sisting of electrically heated tubes through which mineral oil was forced. 
The smoke exited the wing on the under surface very close to the tips at 
about mid-chord. From visual observations, the smoke was essentially 
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completely entrained into the vortex system in a thin annulus at about 
the core radius as long as . the oil flow was steady. If the oil flow 
"chugged", puffs of smoke would appear which would be contained within a 
cylinder of perhaps 5 or 6 core radii concentric with the core, and these 
puffs would be convected downstream without any apparent dissipation or 
radial diffusion within the run length from the test section down to the 
turning vanes of the tunnel. 

The hot-wire probe was mounted on a rotating arm so that it would 
pass through the vortex sufficiently rapidly that the vortex could not 
move significantly during the passage time of the probe. The signals were 
taken out of the arm through slip rings and recorded on magnetic tape. 

The minimum noise level obtained on the tape recorder alone (in a well- 
adjusted condition) with a 2-volt, full-scale range was about 10 millivolts. 
The hot-wire probe was calibrated before each run, since the calibration 
was found to change slightly between runs (from hour to hour) . Calibration 
points at five to ■ ten speeds were taken, since there was some scatter in 
defining the straight line relation between (volts) 2 and (velocity) ^ 2 . 

All calibration and vortex data were taken without smoke so as not to 
foul the wires on the probe. 

Some conditions were experienced which exceeded the acceptance half- 
angle of the probe, which was 35°. The arm was rotated at a speed such 
that the tangential velocity equalled the tunnel air speed. The probe 
was then mounted at 45° to the plane of arm rotation, such that the nominal 
relative velocity was aligned with the probe axis . Under certain combi- 
nations of high wing angle of attack (high T) , probe location and probe 
direction entering the vortex center, the vector sum of the tunnel velo- 
city, the probe translation relative velocity and the peak tangential 
velocity in the vortex produced a velocity making an angle with respect 
to the probe axis as great as 40°. 

During the course of the tests, it was found to be extremely useful 
to have an oscilloscope with a "memory" display feature in parallel with 
the tape recorder and connected to one of the three wires. Since the 
vortex meander in the tunnel (due to turbulence) caused only a small 
proportion of the passes to be core penetrations, the use of the oscillo- 
scope provided an excellent means of determining when core penetrations 
began to occur as the instrumentation tower containing the rotating arm 
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was translated across the tunnel. Later, some oscilloscope photographs 
were taken of the outputs of the wires as recorded on analog tape in order 
to provide data for a manual check of the data reduction program for trans- 
forming the velocity components from the probe to the vortex coordinate 
system. An example of a set of scope photographs is shown in figure 1. 

The data reduction method involved having the analog information on 
the magnetic tapes digitized and entered directly into a computer for 
calculating the corrected velocities . It was found that there is a com- 
bination of minimum tape recorder speed and digitizing rate necessary to 
give adequate definition of the trailing vortex velocity profile. When 
the signals for the vortex, shown in figure 1, were digitized, there were 
seven points in the core. This is probably the minimum allowable definition 
for the tangential velocity profile, but it is inadequate for the axial 
velocity profile. This vortex was generated by the 32-inch span wing, and 
the profile is 55 chords downstream. The core diameter is 0.68 inch. The 
probe traversed through the vortex at 90 feet per second. The tape record 
was generated at 60 inches per second and played back at one-eighth that 
speed. The digitizer cycle time was 1 millisecond. As a check on the 
digitization process , a sine wave was recorded on the tape prior to 
recording the hot-wire data on each run . 

NOMINAL FLIGHT VORTEX CHARACTERISTICS 


Certain decisions have been made by Ames Research Center concerning 
the aircraft to be used and some flight parameters. On the basis of 
these data, calculations have been made to define the wake characteristics. 
The results are summarized in this section. 


The generating aircraft is a T-33 without tip tanks. The pertinent 
characteristics of this aircraft are given in Table I. The flight tests 
are to be conducted at 12,500 feet altitude at a speed of 300 feet per 
second. The results of reference 1 indicate that the T-33 wing planform 
should have a nearly elliptic span load distribution. With this assumption, 
the trailing vortex span is tt/4 of the wing span. The equation relating 
the lift on the bound vortex to the aircraft weight provides the following 
relation for the circulation strength 


r 

o 


4W 

' TTbpU 


(1) 
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For an NACA standard atmosphere (Table II) , the circulation strength of 
the vortices laid down by the T-33 is 819.5 feet 2 per second. The initial 
vortex span is 29.4 feet. 

The specification of the characteristics of a wake vortex follow the 
theoretical model of a vortex developed in reference 2. The nomenclature 
is illustrated in figure 2. The vortex is divided into four regions for 
purposes of specifying the circulation and velocity distribution. There 
is a solid-body-rotation core (the eye) , outside of which are the loga- 
rithmic and outer regions. The outer boundary of the latter is the radius 
(r Q ) at which the vortex strength has achieved 99 percent of its asymptotic 
value. Beyond r Q , the vortex has the classical potential vortex form. 

The initial vortex characteristics after roll-up (some 2 to 4 spans 
downstream of the wing) were estimated on the basis of test results 
obtained in the 40- by 80- and 7- by 10-Foot Wind Tunnels. For the 8-foot 
span rectangular wing at a lift coefficient of about 1, a value of r 2 = 

1.1 inches was measured 55 chords downstream of the wing. The vortex 
radius did not appear to change with distance downstream, as visualized 
by smoke. For tests in the 7- by 10-Foot Wind Tunnel, Corsiglia and 
Chigier found at comparable distances downstream and at the same angle of 
attack that the vortex core radius for a Convair 990 swept-wing model was 
approximately twice that of a rectangular wing. This difference is attri- 
buted to the influence of span load distribution on the rolling -up process, 
since the swept wing had more nearly an elliptical span load distribution 
than the rectangular wing. Further, tests in the 40- by 80-Foot Wind 
Tunnel on the 8-foot span wing and a 32-inch span wing indicated the 
vortex core radius to be approximately proportional to the span. Thus, 
the 1.1-inch radius, scaled up to the T-33 span and doubled to account for 
the influence of span load distribution, yields a value for the flight 
tests of 0.86 foot. This value may be somewhat high but it will be used 
as the nominal T-33 wake vortex radius, r , in the remainder of the 
calculations. y 

As a part of the analysis of data taken from the 40- by 80-Foot Wind 
Tunnel tests, the velocity distributions taken through the core by hot-wire 
probes were examined to obtain the circulation distributions. On the 
basis of these measurements, the radius ratio x l /x was found to be 

about 0.4. Therefore, with T , x^, and r Q known, Tj can be obtained 
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( 2 ) 


from equation (6) of reference 2, 

K I* 

=£ = 0.76 - 0.928 In — 

K . r o 

The following results are obtained: 

T = 819.5 feet 2 per second 
o c 

r = 509 feet 2 per second 

v = 94 feet per second 

r = 0.86 foot 

r =2.15 feet 
o 

These values would obtain 55 chords, or about 370 feet, behind the T-33 
wing. On the basis of the wind-tunnel tests, these values would not 
change appreciably up to 800 feet behind the aircraft. No systematic 
measurements exist for longer distances downstream. 

Axial velocity profiles measured in the 40- by 80-Foot wind Tunnel at 
55 chords downstream of the 32-inch span wing indicate a maximum axial 
velocity defect of about 15 feet per second for a free-stream velocity of 
90 feet per second. The entire defect occurs within a cylinder having a 
radius approximately equal to r Q . Insufficient data have been reduced 
to determine the manner in which this defect changes with free-stream 
velocity and wing lift. It should be noted that the axial velocity profile 
could have large gradients and a complicated structure. The wind-tunnel 
measurements at 55 chords show a profile similar to the sketch 



9 



whereas smoke studies made near the ground show regions near r Q where 
one annulus of flow is moving upstream and one annulus downstream relative 
to an observer on the ground. 

Several things can happen to change the vortex spacing and radius. 

The vortex pair moves vertically downward by self induction. Buoyancy 
effects begin to occur due to altitude changes or engine wake ingestion 
which can have the net effect of continuing the downward motion, causing 
a leveling off of the vortex pair or causing the vortices to rise again, 
depending on the atmospheric temperature profile and the gas temperature 
in the vortices (engine air ingestion). MacCready (ref. 3) indicates 
cases over the ocean where vortices have continued to descend (at least 
at low altitude) whereas most other observations (over land) have indicated 
leveling and/or rising of the vortices after an initial descent. The 
vortex spacing and strength could change the order to 10 to 20 percent due 
to the influence of the buoyancy forces over a typical vortex wake life- 
time (ref. 2). The vortex radius would tend to increase, probably rather 
slowly, with vortex age due to diffusion outward of the vorticity. The 
rate at which this happens depends on atmospheric turbulence and is not 
predictable with confidence at this time. 

The persistence of the vortices also depends greatly on atmospheric 
turbulence. Typical persistency times range from 60 to 150 seconds, which 
in distance correspond to approximately 3-1/2 to 8-1/2 miles at 300 feet 
per second. 

Turbulence tends to promote breakup of the vortex system either by a 
Crow instability, in the event of energy in the turbulence having wave 
lengths comparable to the vortex span, or by single vortex bursting due 
to axial pressure gradients. Persistency of 120 to 140 seconds was 
obtained behind a Boeing 747 over land at 4000 feet altitude (ref. 4) and 
of 40 to 70 seconds behind a DeHavilland DHC-2 at 1000 feet altitude in 
"relatively calm" air over land (ref. 5). In the latter tests, increasing 
turbulence levels changed the dissipation mode from single vortex dissi- 
pation to the Crow instability. Contrail tests reported in reference 6 
on a T-33 aircraft indicate contrail lengths as long as 20,000 feet, al- 
though the speed and altitude of the aircraft are not specified. If 
flights are to be made over the ocean, it is probable that the turbulence 
levels are somewhat lower than over the land, and a reasonable estimate 
for the vortex trail length would be 60 to 90 seconds, or 3-1/2 to 5 miles. 
10 



ATMOSPHERIC TURBULENCE MEASUREMENTS 

It is recognized that atmospheric turbulence plays ' a significant role 
in both the time to dissipation of a vortex and the mechanism of dissipa- 
tion. At present, however, no rational method for predicting the effect 
of turbulence on vortex dissipation has been developed. It is desirable, 
nonetheless, to obtain measurements of atmospheric turbulence levels (as 
well as other atmospheric conditions) during the flight tests to com- 
pletely document the environmental conditions and provide data for subse- 
quent theories which will be able to treat the influence of turbulence. 

Turbulence wave lengths long compared to the generating aircraft 
span can generally be expected to be important in deforming the trailing 
vortices and amplifying instabilities of the Crow type. Very short wave 
length turbulence will determine the turbulent shear, u'v' , at the edge 
of the vortices, and so is important for determining the expansion in 
diameter of each trailing vortex. It is important therefore to determine 
how to measure the atmospheric turbulence, yet keep the effort to be 
expended on its measurement in correct perspective when compared to the 
most important task of obtaining high quality vortex velocity profiles. 

Figure 3 is a power spectrum for atmospheric turbulence. It results 
from a composite of information from references 7, and 8. Only relative 
units are given on the ordinate, as the whole curve will shift up or down 
with increasing or decreasing turbulence intensity. Fortunately, the 
turbulence wave lengths that are of concern for the present study are 
1000 feet and less. 'The turbulent structure in these wave lengths is 
independent of how the much larger eddies were produced, and this range 
is known then as the universal equilibrium range. This range is divided 
into two parts. The upper portion down to 0.1-foot wave length is the 

-5/3 • 

inertial subrange. The energy there decreases as k , where k is 
the inverse of wave length. At shorter wave lengths, where viscous 
dissipation is important, the energy decrease with k is often approxi- 
mated with k 7 (ref. 9). A negligible fraction of the total turbulence 
energy is in this dissipation region, although this turbulence may still 
be important to vortex diffusion. 



Since the shape of the turbulence spectrum in the range of wave 
lengths of interest to the trailing vortex problem is known, but the over- 
all level of the curve is unknown, only one measurement is needed to define 
the curve. The simplest measurement would appear to be the root-mean- 
square of the turbulence signal from one hot wire over a specified fre- 
quency range. The following analysis has been made to determine rms signal 
levels that can be expected and to identify problems with this measurement. 

MacCready (ref. 7) has compared and quantified atmospheric turbulence 
levels with pilot's qualitative reports of its effect on the aircraft, 
i.e., "zero," "light," "moderate," etc. It is generally expected that 
trailing vortices persist longest at the lowest turbulence levels, so the 
turbulence rms levels should be low where the planned flight tests are 
performed. Selecting the limit between "zero" and "light" turbulence from 
Cessna 180 tests, the turbulence spectrum in the inertial subrange (0.1- 
foot to 1000-foot wave lengths) as sensed by an aircraft with velocity U, 
feet per second, can be described by (ref. 7): 

j G(k) = E (k) = 0.585xl0~ 3 U 2/3 f -5/3 (3) 

The frequency is f, Hz; E(k) is the one -dimensional power spectral 
density for longitudinal velocity fluctuations (ft 2 /sec) , and G(k) is 
that for the lateral fluctuations. Although the turbulence is isotropic, 
the velocity fluctuations measured from a moving vehicle are not true 
energies and thus differ slightly for different components (ref. 10) . If 
the aircraft did not respond to this turbulence or the effect of its 
movement were subtracted out, then a single hot wire perpendicular to 
the aircraft axis would measure E(k). 

Integrating equation (3) between two frequencies f x and f 2 gives 
the mean square between those frequencies: 

f 

f E(f)df = - 0.876 U 2/3 (f“ 2/3 - f~ 2/3 ) (4) 

l 

For U = 300 feet per second, f i = 0.2 Hz, and f 2 = 200 Hz, the rms 
turbulence level is 
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As a point of reference, this is a low turbulence level for a wind tunnel. 
Using the method of calculation described immediately below, the turbulence 
level in the 40- by 80-Foot Wind Tunnel has recently been found to be 0.45%. 

Good hot-wire sensitivity and low instrumentation noise is necessary 
to properly measure the above turbulence level. The constant temperature 
hot-wire anemometer sensitivity is determined from its operating equation 

e 2 = A + Bu 1//z (5) 

where e is volts, u is the instantaneous velocity component perpendi- 
cular to the wire, and A and B are operating constants. Larry Olsen 
recently used a DISA hot wire in a test in the Ames 40- by 80-Foot Wind 
Tunnel that is similar to that planned for the flight tests: a 10-micron 

diameter platinum-rhodium wire with plain ends (not gold plated). At an 
overheat ratio of 1.6, he measured 9 volts at 30 feet per second, and 12 
volts at 250 feet per second; thus B =5.57. Differentiating equation (5) 
and rearranging. 



For U = 300 feet per second, e =12.3 volts, and for the above turbulence 
level, the rms signal will be 0.9 millivolts. It would be desirable to 
measure turbulence levels down to about one-third this level of transition 
between "zero" and "light" turbulence, i.e., 0.3 millivolts. 

The DISA 55D35RMS Voltmeter (Large-Scale Aerodynamics has three of 
these) has a much lower frequency limit than other rms meters; its limit 
is 0.1 Hz. Its smallest full-scale reading is 1 millivolt, so the tur- 
bulence level can be measured with this meter. A problem arises, however, 
with the rolloff of an rms meter. For the DISA meter it is conjectured 
to look as presented below. 
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Probable meter 



Hz 

Since the turbulent spectrum energy level keeps rising rapidly below the 
lower frequency limit that is selected, the lower limit should be filtered 
to obtain a characteristic such as shown with the dotted line. Once this 
characteristic is known, the above integration should be redone. There is 
negligible energy above the upper frequency limit, so there is no need 
for a cut-off filter there unless it is to eliminate instrumentation 
noise. 

It should be noted that the noise level of a magnetic tape recorder 
is too great to allow the hot-wire signal to be recorded and the rms value 
determined later. Under good laboratory conditions, the minimum noise 
level for a good tape recorder set for a 2 -volt range (typical of the 
range needed for the above hot-wire sensitivity) will be about 10 milli- 
volts rms. It is therefore recommended that the rms level of at least 
one channel , preferably all three , of the triple sensor hot-wire anemo- 
meter be read directly in the flight test experiment during straight and 
level flight outside of the vortex. (The readings need not be simultaneous, 
but may be switched through one meter . ) The DISA rms meter has an external 
plug for the mean square value, and this output could.be recorded, if 
preferred. 

The noise level on the DISA model 55D05 battery-operated constant 
temperature hot-wire anemometer is advertised to be approximately 0.8 
millivolts rms at full band width. The band width of the anemometer 
depends on the probe wire diameter and measuring conditions, but will 
probably be at least 50 KHz. It is expected that the instrument noise 
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level can be greatly reduced by using a high frequency filter at 5 to 
10 KHz, and thus still make possible the measurement to the above listed 
turbulence levels. . 


The above analysis for the response of the hot-wire anemometer to 
the atmospheric turbulence assumes that the hot wire is mounted on a 
perfect platform that does not respond to the turbulence. The short 
period frequency of the Lear jet is 0.6 Hz, and the boom natural frequency 
is about 11 Hz, both within the integration band width used above. The 
data can be examined after the flight for energy at these frequencies and 
some smoothing can be done if it appears necessary. 

Another possibility that should be considered is operation of the 
data recording system at higher sensitivity levels compatible with meas- 
urement of turbulence velocities for enough periods of time away from the 
vortex wake to attempt to establish ambient turbulence, levels . 


The vanes on the boom also respond to the lateral turbulence compo- 
nents v' and w ' . Again making the simplifying assumption that the plat- 
form on which the vanes are mounted does not respond to the turbulence, as 
for the previous analysis, the output signal from the vanes is easily 
determined. Denoting the perturbation angle of the flow due to turbulence 
in the vertical direction as a', then 



tan 


w' 

U + u' 



(7) 


The vanes only respond accurately to the frequencies up to 15 Hz. Changing, 
the upper frequency limits in the integration performed in connection with 
equation (4) from 200 to 15 Hz reduces the integral by 2 percent. Including 
the 4/3 factor for lateral turbulence. 



0.0014 radians 


For the vane angle transducer sensitivity of one volt per 10°, there 
results 



millivolt rms 
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Dur.ing the series of tests on trailing vortices in the 40- by 8 0- and 
the 7- by 10-Foot Wind Tunnels, it was observed that the vortices were 
undergoing significant lateral oscillations which appeared to increase in 
amplitude with increasing distance downstream of the wing. This behavior 
was believed to be induced by turbulence in the tunnels. Motion pictures 
were taken of the vortex motion and some statistical analysis was done on 
the photography. In addition, a theory was developed to explain the 
results. The details of this work are reported in reference 11. Since 
the work yields some indication of the anticipated random motion of the 
vortices in flight, the major results are reported here. 

The motion of the vortices in the tunnels was found to be quite 
random. The correlation between the vertical and horizontal deviations 
from the mean position of a vortex was found to be quite high, however, 
which might indicate that some organized form of motion such as a "cork- 
screw" form was occurring. The rms values of the displacement of the 
vortex from its mean position (df 2 ") were calculated at several locations 
and are shown as the data points in figures 4(a) and 4(b), for the 40- by 
80-and 7- by 10-Foot Wind Tunnel tests respectively. The turbulence 
velocity levels, a, are measured values for the two tunnels. The data 
in figure 4 (a) are for both the 8-foot and the 32-inch span wings and show 
that the rms displacement is independent of the wing scale. The reason 
for the difference between right and left vortices is not known, but may 
be a result of a consistent variation across the tunnel of turbulence 
level. 

A theory was developed for the vortex motion at a distance d down- 
stream of the wing, based on the premises that the vortex is a filament 
which has the transverse velocity locally of the turbulence at that point 
and the vortex does not respond to high frequency turbulence. The results 
of the analysis indicate that the root-mean-square deviation of the vortex 
varies with d as shown. 



where A^ is a scale wave length of the turbulence. For cases where 
A /d is large (small d or large A ), the arctan is tt/2 and 

S S 
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Equation (8) is represented as the lines in figure 4, with the values of 
the parameters shown. The values of the turbulence scale selected were 
values approximating the cross section dimensions of the tunnels, which 
are shown to cause reasonably good agreement with the data. 

The theory can be applied to the flight test conditions, provided 
the form of the power spectrum is known. Using the model from reference 7 
for the inertial subrange, the results of figure 5 were obtained. These 
results indicate that the rms displacements reach the order of the wing 
semispan for light turbulence with typical values of the turbulence length 
scale at distances of 2 to 4 miles aft of the aircraft. Provided that 
single vortex bursting has not occurred within this distance, this calcu- 
lation would indicate the minimum distance required to permit the 
"touching" of the two vortices and the subsequent formation of vortex 
rings, in accordance with the Crow instability model. 

DATA REQUIREMENTS FOR THE FLIGHT TESTS 

A principal objective of the flight test program is to obtain velocity 
profile measurements of wake vortices as an input to the problem of calcu- 
lating the severity of aircraft-wake encounters. In addition, it is 
desirable to obtain sufficiently detailed measurements of both the vortices 
and the environment to permit evaluation of theories describing the for- 
mation, motion, and dissipation of wake vortices. For this purpose, 
existing theories have been examined to determine the kinds of data 
required either as input or as a means of evaluating the predictive cap- 
ability. The theories examined were the following: 

(a) The NEAR, Inc. vortex structure theory, reference 2. 

(b) The work of McCormick, reference 12. 

(c) The NASA work of Baldwin and Chigier, reference 13. 

(d) The instability theory of Crow, reference 14. 

(e) The NEAR, Inc. theory on buoyancy effects, reference 2. 
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(f) The buoyancy work of Scorer and Davenport, reference 15. 

(g) The. buoyancy work of Tombach, reference 16. 

The data fall into three categories : airplane data , wake vortex data 

and meteorological data. The description of these data, together with a 
discussion of means of obtaining .the data, is given below. 

Airplane Data 

The airplane data required are those necessary to define the circu- 
lation strength of the wake vortex system of the generating aircraft (the 
T-33) . In accordance with equation (1) , the following information is 
required : 

(1) Aircraft weight at the time the wake is being generated. 

(2) Aircraft air speed during a probe period. 

(3) Air density at the altitude of the generating aircraft. 

(4) Some indication of the constancy of the flight path of the 

generating aircraft. The ideal quantity would be a recording 
of vertical acceleration. The pilot's assessment would pro- 
bably be satisfactory. For example, it has been shown (ref. 5) 
that small pitch oscillations of an aircraft at the correct 
frequency can decrease the time to onset of the Crow instability 
by a factor of 2. The object in the ARC tests is to maintain a 
constant wing lift (e.g., no vertical accelerations) so as to 
lay down a uniform wake. 

(5) Some indication of the accuracy with which the air speed is 

maintained during a run. If the quantity is recorded, the 
records would provide an indication of the excursions. Again, 
the pilot's assessment would probably be satisfactory since 
this is not a critical area. 

Trailing Vortex Data 

The data required on the trailing vortices consist of velocity and 
temperature traverses through the vortices, auxiliary information to fix 
the location of the traversing station relative to the generating aircraft, 
and some photographic coverage. The considerations involved are discussed 
below. 
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The principal measurements to be obtained are three components of 
velocity through the vortex. These are to be obtained with a three-wire 
DISA hot-wire probe mounted on the nose boom of the Lear jet. From the 
NEAR, Inc. hot-wire experience in the 40- by 80-Foot Wind Tunnel tests, 
several things are worthy of comment. 

A portion of the core flow in the vortices may be made up of engine 
exhaust air, which is hotter than the ambient and therefore a potential 
source of error for the hot-wire measurements. The rate at which the 
engine exhaust mixes with the ambient air in the presence of the vortex 
system is not predictable at this time. A rough estimate was made 
assuming only a coflowing turbulent jet exists, and the calculations 
indicated that a 170 to 1 dilution of exhaust gases would occur within 
approximately 1700 engine exhaust radii, which is about one mile. If 
the engine exhaust temperature is 1000° F at the tailpipe, the mixed gas 
temperature is about 20° F at one mile, or some 6° above ambient. Thus 
a reasonable part of the vortex wake region of interest may have tempera- 
ture differences above ambient that must be accounted for in the hot-wire 
data reduction, since 1° C temperature difference causes an error in 
velocity of about 2 percent. It is also of interest to know what the 
temperature profiles are in the vortex in order to apply buoyancy theories 
to the vertical motion of the vortex pair. 

With respect to correction of the hot-wire results for temperature, 
Bradshaw (ref. 17) recommends that the calibrations be plotted in terms 
of Nusselt number (Nu) versus Reynolds number (Re) for the temperature 
conditions of the calibration procedure. Operation at other temperatures 
can then be evaluated through use of theoretical behavior for the effect 
of temperature on these parameters. The Reynolds number of the flow over 
a 10-micron wire at standard atmospheric conditions at 12,500 feet reaches 
a maximum of 45 in the vortex. This value is just at the transition 
region where the wake behind a cylinder begins to oscillate and affect 
the heat transfer. For a Reynolds number up to 44, Bradshaw recommends 
the heat transfer relation 



17 

= 0.24 + 0.56 Re 0 - 45 


( 10 ) 


19 


where 


Nu = H/iT£k(T - T ) 

' w a 

and 

H = heat transfer rate to the fluid 

£ = wire length 

k = conductivity of air 

T w = wire temperature 

T = air temperature 

cl 

T m = average of wire and air temperature 
and Re and Nu are evaluated at T . 

The Lear jet carries a single wire probe for the measurment of tem- 
perature. Since this wire must have sufficient response to follow vari- 
ations through the vortex, its thermal inertia characteristics were 
calculated to determine if a d.c. bridge might offer satisfactory response 
Using the same wire type as for the three-wire probe, the wire time con- 
stant is about 3.5 milliseconds. However, for a 5-percent error in 
following a sinusoidal temperature oscillation, the maximum allowable 
frequency is only 15 Hz, which is too low for the vortex penetration case. 
Therefore, constant-current hot-wire anemometer electronics will be neces- 
sary to amplify and compensate the wire signal to obtain reasonable 
temperature frequency response (1000 Hz is adequate) . 

The acceptance angle of the hot-wire probe was a problem in the 
wind-tunnel tests. In the flight tests, however, with a Lear jet velocity 
of 300 feet per second and a maximum tangential velocity in the vortex 
of 94 feet per second, the flow angles on the probe should not exceed 
17.4°, which is well below the acceptance angle of 35°. 

It should be noted that traversing directly through the center of an 
axisymmetric vortex will not produce a symmetric trace on a wire of the 
three-wire probe, even if it is aligned in the plane perpendicular to 
the vortex. For instance, for the above peak tangential velocity and 
traverse speed, and for typical values of A = 50 and B = 5.5 (eq. (5)), 
the voltage of the larger peak (or deflection on an oscilloscope) should 
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be 1.79 times that of the smaller peak. This is due to two factors that 
add, the nonlinearity between velocity and voltage in a hot-wire charac- 
teristic, and the nonlinear effect of angle change on the hot wire, as 
sketched here. The velocity components perpendicular to the wire for 
0 = 17.4°, 0°, and -17.4° for U = 300 feet per second are 232,245 and 
300 feet per second, respectively. 

0 = 17.4° 



These velocities with the nonlinear velocity-voltage characteristic 
provide the unequal peak voltage measurements . 

Calibration of the probe should be done at five to ten velocities in 
order to establish the voltage-velocity calibration curve, since there is 
always some scatter in the data. It would also be desirable to calibrate 
several times during a flight because of the possibility of smoke or 
particle impact changing the wire characteristics. 

It was noted earlier that an oscilloscope with a "memory" feature 
was used to monitor the output of the wire that most nearly sensed tan- 
gential velocity in the wind-tunnel tests. It is recommended that this 
be done in the Learjet in order to determine when core penetrations are 
made, since so much more information is available with a core penetration 
than without. 

Further, it was noted earlier that high tape recording speeds were 
used in the tunnel tests so that the tape could be slowed in the digi- 
tizing process to gain definition in the velocity profile in the core 
region. For flight test conditions of 300 feet per second, a core dia- 
meter of 1.7 feet, a digitizing cycle time of 2 milliseconds and no 
tape speed reduction or playback, the number of data points in the core 
would be only 2.9. Therefore, a 2 to 1 tape slowdown is absolutely 
necessary, and a 4 to 1 slowdown is desirable, particularly for measuring 
the axial velocity profile. If the minimum digitizing speed is 1-7/8 
ips, then a recording speed of 7-1/2 ips is preferable when flying 
through the vortices. 



Another possible limitation on the tape recording is the frequency 
response. The basic sine wave component of the tangential velocity pro- 
file of the trailing vortex for the above conditions is 87 Hz. For 1- 
percent accuracy in recording this signal, the recorder must be capable 
of recording 13 times this fundamental frequency, or 1.1 KHz. An fm tape 
recorder operating with a carrier frequency of 13.5 KHz (standard) at 
7-1/2 ips has a frequency limit of 1.25 KHz. Thus 7-1/2 ips appears to 
be the minimum acceptable tape speed for recording vortex information 
unless the digitizing cycle period is reduced and the analog tape 
recorder has extended band width. 

It is necessary to determine the location of the wake traverse rela- 
tive to the generating aircraft. It is understood that some airborne DME 
equipment can be obtained from Langley Research Center which will provide 
a separation distance measurement to within about 30 feet. It is desirable 
also to measure the vertical distance that the vortices have moved between 
the time they were laid down and the time the Lear jet passes through the 
vortices. The objective here is to obtain information on the vertical 
translation of the vortices to evaluate buoyancy theories. The vertical 
motion of the vortices should be the order of several hundred to one. 
thousand feet maximum. If the T-33 is able to fly at a perfectly constant 
altitude, the vertical translation is the difference in altitude between 
the two aircraft at the time of vortex entry by the Lear jet. Since each 
altitude is measured individually and differenced, measurement of each 
altitude to +10 feet will provide +20 feet on the difference or about 5- 
to 10-percent accuracy which should be sufficient. Since the difference 
is the major feature of interest, the altimeters in the two aircraft can 
be adjusted to read identically with the aircraft side-by-side prior to 
beginning testing. In order to ascertain whether the T-33 is flying a 
constant altitude, the pilot can radio altimeter readings to the Lear jet 
for voice recording periodically, so that his altitude history can be 
reconstructed. 

With regard to the number and location of the axial stations which 
are probed, it would appear desirable to have one station as close to the 
generating aircraft as is feasible (perhaps 1/2 mile.) ,. three or four 
stations in the stable portion of the vortex, and one or more in the 
breakup region. This would amount to spacings between probe stations of 
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the order of one-half to one mile depending on the distance to vortex 
breakup. The last station (breakup) is deserving of some comment. 

From the standpoint of alleviating the wake vortex problem, an under- 
standing of the dissipation mechanism is of great importance. The 
instability mechanism predicted by Crow (ref. 14) agrees with observations 
of vortex breakup, but his theory does not postulate how the instability 
develops, what influence atmospheric turbulence has on the development, 
nor what happens when the flow gets highly nonlinear in the breakup and 
ring formation region. No measurements have been made of the velocity 
distributions in the ring formation region to determine the potential 
hazard to passage of aircraft. Also, in the case of single vortex dissi- 
pation, no measurements have been made to indicate whether a bursting or, 
a more rapid vorticity diffusion process is responsible for the 
dissipation. 

We would suggest that consideration be given to making some traverses 
in the dissipation region. The exact approach to how these should be made 
would probably have to be decided upon during the test program when the 
dissipation mechanism has been observed. It would probably be desirable 
to make traverses just upstream of the dissipation event and within the 
region where dissipation has occurred. 

Certain photographic coverage is considered essential. The first 
category is coverage to determine the development of the Crow instability 
and the dissipation mechanism. The wave length of the vortex oscillations 
for maximum growth rate of the Crow instability is 6.6 spans, or about 
200 feet. It is recommended that the Learjet or another chase plane be 
flown parallel to the vortex system at 300 feet per second either just 
above or. just below the vortex trail and movies be taken looking forward 
and slightly down (or up) . The vortex lateral motion as derived from 
these movies will represent the time variation of the vortex oscillations 
at a fixed distance behind the aircraft. It is suggested that this 
coverage be obtained at three or four axial stations behind the generating 
aircraft for cases where the dissipation is due to a Crow instability and 
where it. is due to single vortex bursting or diffusion. A single camera 
can be used, and the variation of the separation distance between the 
vortices will appear very close to its true perspective in the photography. 
Measurements from the photographs can be related to a time average of this 
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distance to obtain fractional deviations and their frequency. It is 
recognized that this photography may be difficult to obtain due to overall 
deformation or twisting of the vortex system due to large scale turbulence, 
buoyancy effects, etc. Nevertheless, little, if any, such systematic data 
has been obtained, and it is worth some effort to collect such data. 

Other categories of photographic coverage would include photography 
to document the form and details of vortex dissipation and the influence 
of the Learjet on the vortex system at and shortly after the instant of 
passage of the aircraft through the vortex. The former is of interest in 
order to gain greater understanding of the details of the dissipation 
region. Coverage would be desirable both from a side view and a top view. 
The coverage of the Learjet penetration is of interest in terms of asses- 
sing the influence of the aircraft on the stability of the vortex. On 
the basis of the 40- by 80-Foot Wind Tunnel tests, it is expected that the 
Learjet should not cause vortex breakup, since the size of the "probe" 
relative to the vortex core is not greatly different in the two cases. 

Meteorological Data 

It is essential to document the meteorological conditions in which 
the vortex measurements are made, since they have an important part in 
determining the motion and dissipation of the vortex system. The consi- 
derations involved in measuring turbulence were discussed in a previous 
section. In addition to these data, information is desirable on the 
temperature and pressure at the test altitude, and the manner in which 
the temperature varies with altitude within, say, 500 feet above to 1500 
feet below the test altitude; i.e., the range over which the vortex system 
might be expected to translate vertically. 

From the standpoint of calculating buoyancy effects on the vortex 
vertical translation, it would be desirable to know whether or not vertical 
velocities are present in the atmosphere. The vertical velocity components 
with wave lengths of 1000 feet or less can be obtained from the turbulence 
measurements, with the assumption of isotropy, which is probably a reason- 
able assumption. The longer wave lengths, which in the limit become a 
uniform vertical velocity (at least over the length scale of the vortex) , 
are difficult to measure because the aircraft tends to follow these. In 
addition, since the mutual induction velocity of the vortex pair is the 
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order of 5 feet per second, the long wave length vertical velocities of 
interest are small compared to the forward speed of the aircraft. It may 
be possible by examining angle of attack, altitude, and vertical accele- 
ration data outside the vortex to back out long wave length vertical air 
velocity data, and it is suggested this be considered in connection with 
any buoyancy calculations made with the data. 

Another measurement of considerable potential importance is humidity. 
In a subsequent section of this report, the influence of condensation on 
vortex behavior is discussed. It is noted there that tests at low alti- 
tudes under humid and/or hazy conditions (ref. 5) have shown distances to 
vortex breakup of one-third those obtained in clear weather. While these 
tests did not identify the mechanism responsible for the early breakup, 
it is possible that high humidity is responsible, at least in part. There- 
fore it is recommended that humidity measurements be made on board the 
Learjet at several times throughout a flight test and within the vortex, 
if possible. The latter location is of interest because of the possibility 
of high water vapor content due to ingestion of engine exhaust air. 

INTERACTION BETWEEN THE PROBE AIRCRAFT AND THE VORTEX WAKE 

The flow disturbance around the aircraft in steady flight will induce 
velocity components at the probe which will appear as biases changing with 
angle of attack. Also, the probe aircraft will respond to the flow field 
induced by the vortex wake, and the resulting motion will induce a velo- 
city relative to the hot-wire probe which must be removed from the data 
and could cause problems in aiming the aircraft so as to penetrate the 
core. Finally, the presence of the aircraft will cause the vortex to move 
or deform, which could affect the measurements. Some estimates of the 
sizes of these effects were made to attempt to anticipate problems in 
this area. 

Steady Flow Disturbance Due to the Aircraft 

The perturbation flow field around the Learjet in a steady flow at 
several angles of attack was calculated using an existing potential flow 
program for a wing-body combination. The program is limited to an axi- 
symmetric body shape. Consequently, the correct area distribution was 
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used, but the Lear jet fuselage was modeled as an axisymmetric body. The 
program fits a source and doublet distribution to the fuselage and a vortex 
lattice to the wing. The tip tanks were not modeled because their influ- 
ence should be small. The details of the configurations were obtained 
from reference 18. 

Results for the perturbation velocities over the forward part of the 
Lear jet are shown in figures 6 and 7. Figure 6 shows longitudinal and 
vertical velocity components in the vertical plane of symmetry for 0° and 
10° angle of attack at various distances above and below the fuselage 
centerline. The vertical distance scale has been greatly expanded to 
simplify the presentation. At the hot-wire location (x = 5 feet) , the 
axial (u ) component is about 3 percent of free stream (6 feet per second) 
and essentially does not change with angle of attack, since it is due to 
a fuselage-thickness effect. The vertical component ( w a ) is zero at 
a = 0 and increases to about 1/2 percent of free stream (1.5 feet per 
second) at a = 10 degrees due to a wing-induced upwash. These velocities 
are small compared to the peak tangential velocities in the vortex. In 
particular, if the angle of attack is perturbed due to vortex-induced 
motion of the aircraft, the only influenced component will be the upwash 
velocity, whose magnitude changes so little that - its effect on the probe 
readings can be neglected. Figure 7 shows u and w velocities in a 
vertical plane 10-1/4 feet from the plane of symmetry (60-percent span 
station) . These results are of interest principally in estimating vortex 
deformation due to the passage of the aircraft and show the upwash forward 
of, and downwash behind, the wing and the acceleration (- u ) over the 
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wing and deceleration under the wing. 

Aircraft Motion Induced by the Vortex System 

The Learjet will respond to the flow field induced by the vortex 
system. The potential problems are control of the attitude of the air- 
craft and the relative velocities induced at the hot-wire probe by air- 
craft motion. For the case where the aircraft flies along the vortex and 
penetrates by descending or ascending through the vortex, the worst case 
is that in which the aircraft finds itself in the center of one vortex 
with a large induced roll, as shown in the sketch. 
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An estimate of the induced rolling moment due to both vortices acting on 
the aircraft was made assuming a wing alone to be centered on one vortex. 
Potential vortices with solid body rotation inside r x were used to 
obtain spanwise upwash distributions along the wing panels, and strip 
theory was used with a uniform section lift coefficient along the span to 
obtain the vortex- induced lift distribution and the resultant rolling 
moment. A value of = 0.3 was calculated. The control roll moment 
available through differential aileron deflection was estimated assuming 
20° of aileron deflection, a change in section angle of zero lift with 
change in aileron deflection angle of 0.6, and a section lift coefficient 
slope of 0.07 per degree (equal to the wing lift coefficient slope). 

The resulting rolling moment coefficient, calculated again using strip 
theory, is = 0.36. If the assumptions are reasonable, these results 
indicate the Lear jet does have sufficient control effectiveness to main- 
tain a zero bank angle in the vortex, but the two values are so close that 
there are likely to be large excursions in roll attitude as the aircraft 
enters the vortex in spite of the skill of the pilot.. Although the atti- 
tude of the aircraft will be monitored so as to extract the aircraft 
dynamics from the hot-wire signals, it is felt that the marginal control 
situation is not conducive to satisfactorily controlling the trajectory 
of the aircraft as it makes the vortex penetration. 

For a case in which the aircraft approaches the vortex system hori- 
zontally at right angles to the wake, the aircraft will pitch and heave 
as it passes through the vortex. An estimate of the vertical acceleration 
was made assuming the wake-vortex-induced lift on the Lear jet to be given 
by the vertical velocity induced at the wing quarter-mean-aerodynamic- 
chord by the two vortices ,' combining this velocity with U to get an 
angle of attack (see fig. 8(a)), and causing the resulting perturbation 
wing lift to act on the aircraft weight. A plot of the resulting accele- 
ration is shown in figure 8 (b) . The values are low and should not cause any 
problems on aircraft loading. The acceleration was integrated to obtain 
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vertical- velocities in order to estimate the difficulty that the pilot 
might have in adjusting his height to hit the center of the vortex. The 
Lear jet would sustain an upward velocity component of only 0.65 foot per 
second by the time the probe reached the vortex. Since the accelerations 
are significant only within 50 feet of the vortex, the aircraft height 
relative to the vortex should change, for a speed of 300 feet per second, 
only a few inches at most due to the vortex-induced upwash. 

An estimate was also made of the pitching rate induced on the aircraft 
by its interaction with the vortex. Two approaches were examined. In 
the first, it was assumed that a pitching moment was imposed equal to the 
moment on the Learjet in a uniform flow at an angle of attack equal to the 
local angle at the quarter-mean-aerodynamic-chord (fig. 8(a)). From ref- 
erence 18, the pitching moment coefficient slope of the complete aircraft 
about c/4 is -0.02 per degree. For the center of mass at the same loca- 
tion and a moment of inertia in pitch of 20,000 pounds-feet per second 2 , an 
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angular velocity of -9x10 degrees per second (nose down) is reached at 
the time the nose boom intersects the vortex. This value is highly depen- 
dent on the assumed center of mass but is so small that any effects due 
to rotation are negligible. The other approach taken assumed that the 
only region of the aircraft acted upon by the vortex-induced upwash is 
the nose of the aircraft, since the gradient in upwash is high near the 
vortex and drops off rapidly (fig. 8(a)). Slender-body theory gives the 
lift curve slope on an expanding body as 2, based on the maximum cross 
section area. Applying this concept to the Learjet with the assumption 
that this lift acts 15 feet forward of c/4, one computes an angular velo- 
city of 0.5 degrees per second when the boom probe enters the vortex core. 
This value is also negligibly small in terms of either causing the aircraft 
to miss the core during the vortex penetration or causing relative velo- 
cities at the probe. 

Effect of the Aircraft on Vortex Deformation and Stability 

The penetration of the Learjet into the vortex wake of the T-33 has 
one and possibly two effects on the vortices: deformation of the line 

vortices and vortex bursting. No calculations of these phenomena were 
considered justified because of their difficulty. However, some qualita- 
tive observations can be made. 
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In the case of a horizontal probe path normal to the vortices, the 
results of figure 6 indicate that there is insufficient disturbance ahead 
of the Learjet to cause any vortex motion or deformation prior to the 
entry of the hot-wire probe into the near vortex. Assuming that both 
vortices still lie in a horizontal plane and measurements through both 
vortices are possible during the passage of the aircraft, the near vortex 
will not be affected appreciably until the hot-wire probe is some 3 or 4 
feet past the near vortex center. The near vortex will then deform to 
pass over, or possibly under, the aircraft in accordance with the pertur- 
bation velocity field induced by the aircraft. The deformed vortex will 
in turn induce some deformation of the far vortex before the hot-wire 
probe enters the far vortex. Since the vortex span is 29 feet, the near 
vortex will be in the region of the wing trailing edge when the probe 
enters the far vortex and will have been displaced about 3 feet vertically. 
It is felt that this deformation should not cause a large effect on the 
far vortex, and the chances for obtaining good data on the far vortex are 
reasonably good. 

Insofar as bursting is concerned, the tests in the 40- by 80-Foot 
Wind Tunnel indicated that passage of the probe through the vortex caused 
no observable oscillations or bursting of the vortex. The probe in that 
case had a diameter of about half the core diameter, whereas in the flight 
test case, the "probe" (the Learjet) is about four times the core diameter. 
It is possible that the vortex deformation will set up disturbances which 
would enhance the onset of the Crow instability after passage of the 
aircraft. 

For a coaxial probe in which the aircraft ascends or descends through 
one vortex, the influence of the aircraft on the vortex at the hot-wire 
location is more difficult' to predict. Although the aircraft disturbance 
at the hot-wire location is very small, it is possible that the presence 
of the fuselage downstream could cause a large disturbance in the vortex 
which could propagate upstream to affect the flow at the hot-wire location. 
We are not aware of any data which would shed any light on this problem. 
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FLIGHT TEST PROBE PROCEDURE 

The two methods of probing the vortex wake that are under conside- 
ration are a horizontal traverse at right angle to the vortex, as illus- 
trated on the left iselow, and a vertical traverse down (or up) 




through the vortex as shown on the right. In the first case the Lear jet 
must fly faster than the T-33 on the upwind legs so that when it turns to 
cross the vortex wake, the intersection with the wake occurs at about the 
desired distance aft of the T-33 on each pass. In the second case, the 
Learjet may fly at essentially a constant speed just over 300 feet per 
second to keep the intersection point a given distance behind the T-33. 

The navigation problem is more difficult with the crossing technique. 

It would appear that the crossing procedure is the more desirable of 
the two for three reasons. First, it appears from the controllability 
calculation previously described that it may be difficult to control the 
attitude and trajectory of the Learjet if it has to fly into the core 
aligned with it because of the large induced rolling moment. Secondly, 
it is probhbly easier to obtain a core penetration with a crossing tra- 
jectory because the vortex will appear relatively stationary to the pilot, 
whereas with the aligned probe, the vortex will appear to move from side 
to side beneath the aircraft due to standing waves in the vortex. Thirdly 
a crossing probe will provide some information on both vortices in one 
pass, whereas the aligned probe will provide information on only one 
vortex. These advantages would appear to outweigh the disadvantages in 
the crossing probe of the more difficult navigation problem and a greater 
time between probes. It is suggested, however, that some of the probes 
in the vortex breakup region be obtained with the aligned probe technique. 
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or possibly the Lear jet be flow at a constant altitude in the breakup 
region, since this procedure would provide some qualitative pilot reaction 
to the degree of turbulence in this region as well as quantitative data. 

The number of axial stations along the wake where data are to be 
obtained was discussed in a previous section. It is suggested that suffi- 
cient passes be made at a given nominal distance behind the T-33 to obtain 
at least one core penetration and, better two penetrations for each dis- 
tance. The oscilloscope records would provide indications of a core 
penetration. 

EFFECT OF CONDENSATION ON VORTEX CHARACTERISTICS 

Condensation trails (contrails) have long been recognized as one of 
the most vivid means of visualizing trailing vortex wakes from aircraft at 
high altitudes. A considerable amount of research was conducted in the 
1950's and early 1960 ' s on means of suppressing contrails (refs. 6 and 19, 
for instance) . This work was concerned principally with the physics of 
condensation and sublimation rather than the fluid mechanics of vortices, 
and consequently has shed little light on the influence of condensation 
on the structure and persistence of vortex wakes. Recently, however, 
Chevalier (ref. 5) reported on some low altitude flight wake vortex experi- 
ments in which he noted on days of high haze and/or moisture content that 
vortex breakup occurred in about one-third the distance required on clear 
days. While no measurements were made of the meteorological conditions 
(temperature, humidity, turbulence , etc .) , the dominant effect seemed to 
be humidity. At the present time there is no theory, or even a qualita- 
tive model, to explain this effect. 

It is well known from Hilsch tube experiments that an energy sepa- 
ration occurs in rotating flows. The fluid in the center has a lower 
static pressure and total temperature than that in the outer flow. If 
moisture is present, the possible effects of moisture are a 'condensation 
in the center of the vortex due to the drop in static - temperature, a heat 
release due to condensation, transfer of momentum from the gas to the 
liquid droplets, and a radial flow of this momentum in the liquid particles. 
Some preliminary calculations were made of these effects in order to com- 
pare the T-33 flight conditions with other flight conditions where certain 



observations have been made and to determine if some mechanism of insta- 
bility could be found. Three cases were compared: the T-33 at 12,500 

feet, a DCH-2 aircraft at 1,000 feet (Chevalier's experiment, ref. 5) and 
a B-47 at 35,000 feet (experiments reported in ref. 6 demonstrating con- 
trail lengths of up to 120 seconds) . The pertinent conditions are listed 
in Table III. These calculations were made ignoring the heat and moisture 
input from the jet engine exhaust (which was not present in Chevalier's 
experiment) which is too difficult to include at the present time. 

An excellent summary of the equations governing the flow in a vortex 
is given in reference 20. In order to make the calculations tractable, 
the form of these equations for negligible temperature variations and 
compressibility effects and for zero fluctuations (laminar flow) was used 
(eq. (3) of ref. 20) . These equations would apply to the core but not 
necessarily to the outer flow which is turbulent. 

The radial static pressure gradient in the vortex is given by 


3P = 
3r 



( 11 ) 


Expressions for the variation of the velocity with radius are given in ref- 

\ 

erence 2. To minimize the computation time, the vortex was divided into 
three regions: a core with outer radius r^, a portion of the logarithmic 

region from r^ to r Jf and a potential vortex outside of r^ Expres- 
sions for the velocity variation in these regions are available in ref- 
erence 2, as follows 


0 < r < r^ 



1.47 





r . < r 
i 


< r. 


* 7 = i + °- 928 ln (^) > 


( 12 ) 


r t < r < «= 



J 


Using these expressions, the static pressure difference between the vortex 
center and the free stream can be calculated from equation (11) . 
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The static temperature distribution in the vortex was estimated using 
the energy equation. This has the form (ref. 20) 
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where I is the total enthalpy, 

1 q = c p T + j (u 1 2 + v 2 + w 2 ) (14) 

In order to use this expression, the following assumptions were made: 

(1) the flow is steady; i.e., 3/3t = 0 

(2) the Prandtl number is one 

(3) the terms involving u (the radial velocity) and 3/3z (axial 

gradients) are small and can be neglected. 

These assumptions reduce equation (13) to the following 


r ^ - v 2 = 0 


Again through the use of equation (12) , the variation of total enthalpy 
through the vortex can be obtained. One interesting result of these cal- 
culations is that the variation in total enthalpy in the core region is 
due entirely to the velocity variation, and the static temperature is 
constant . 

Condensation can occur within the vortex if the relative humidity of 
the ambient air in which the vortex is formed is sufficiently high that 
the drop in temperature in the vortex would tend to increase the relative 
humidity above one. The absolute humidity co is 


m + m 
w a 


1 + (P-./PJ < M =/ M J 

aw a W 
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where p a and p w are the partial pressures of air and water vapor, 
respectively, in the mixture. The absolute humidity at saturation is the 
value of to for p = p . (T) . Thus, if the ambient air is saturated at 
p . (T ) and T, is lower than T , tu , (T,) will be lower than go . (T ) 
and the difference represents the amount of water vapor which must be con- 
densed when the (saturated) ambient air is entrained in the core as the 
vortex is formed. Saturation pressures at low temperatures were obtained 
from reference 15 and are reproduced in figure 9. Since the predicted 
temperature differences between ambient and the core are small, the slopes 
of the saturation pressure curve of figure 9 were used to calculate changes 
in humidity at the nominal temperatures of interest. 

If condensation does occur within the vortex, heat is released and 
tends to warm the air-water vapor mixture in the core. When the amount of 
water vapor condensed per pound of mixture is known, the heat of evaporation 
and the specific heat of the gas mixture can be used to calculate a tempe- 
rature rise in the core gas mixture. 

If water vapor condenses in the core and is centrifuged outward due 
to its higher density, the core suffers a loss in momentum. The amount of 
the loss depends on the velocity imparted to the condensed water. The 
water droplets will assume tangential and axial velocities equal (essen- 
tially) to the values of the air at the point where the vapor condenses. 

If the centrifugal separation occurs slowly enough, the droplets will 
assume the local tangential velocity as they move outward. Assuming this 
to be the case, the angular momentum per unit mass of core flow assumed 
by the droplets is the product of their mass per mass of core flow and 
v . One can compute the total angular momentum of the air in the core per 
mass of air in the core as 

f 1 p.uv dr 

core ang. mom. _ o ■ . 

core mass r t 

f pu dr 
o 

For purpose of making this calculation, u was considered constant within 
r x . , The ratio of the droplet momentum to the core momentum can then be 
obtained. Using equations (12) for the velocity distribution in the core 
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and logarithmic region out to r , one obtains for the momentum ratio 
1.53 times the mass of vapor condensed per mass of core mixture. 

The results of the calculations described above are shown in Table IV 
for the three conditions noted previously. The DCH-2 has a low wing 
loading compared to the other two aircraft and therefore a low vortex 
strength (Table III) . Even though the ambient pressure and temperature 
are highest for this case because of the low altitude, the low vortex 
strength causes only small changes in pressure and temperature between the 
vortex center and ambient. As a result, the various calculated effects of 
condensation are very small compared to the other two cases. 

The B-47E has the highest vortex strength and therefore the greatest 
relative influence on condensation in the core. However, there is an 
order of magnitude less water vapor in the saturated condition at 35,000 
feet than at 12,500 feet. Consequently, the absolute magnitude of the 
condensation effects is smaller than for the T-33 case. It is of interest 
to note that Pilie (ref. 19) gives a lower limit of water droplet concen- 
tration for contrail visibility of 10~ 5 pounds per pound. While this 
figure may be applicable only under ideal conditions and therefore some- 
what low, the three cases show condensed water concentrations of this order. 
Finally, it should be noted that in the B-47E case, the ambient water vapor 
level is so low that the water content in the jet engine exhaust, which 
was neglected in these calculations, is going to have a very large effect 
on condensation in the wake vortex system. 

These calculated results are not conclusive. In the DCH-2 case where 
Chevalier (ref. 5) observed a much earlier vortex dissipation in moist 
conditions (20 seconds) than in dry conditions (1 to 2 minutes) and where 
the engine (a piston engine) has little moisture input to the vortex wake, 
the calculated effects of condensation are so small that no obvious 
mechanism of dissipation is apparent. In the B-47E case where sufficient 
condensation has been observed to occur to cause the wake to be visible 
for up to 15 miles or so, condensation does not appear to be a dominant 
factor in vortex dissipation. Since Chevalier did not make any measure- 
ments of either the ambient atmospheric conditions or the vortex wake, 
it is not possible to determine the precise mechanism acting on his 
vortices under his "early morning fog or haze" condition which caused early 
dissipation. It is felt, however, that the T-33 case will be more like 
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the B-47E case than Chevalier's case for several reasons: 

• The tests to be performed with the T-33 will be made under 

conditions of good visibility (i.e., low humidity) for safety 
purposes and for purposes of locating the wake so the Learjet 
can effect wake penetrations ^ 

• The work reported in reference 6 discussed contrail lengths 

behind T-33 aircraft of up to 20,000 feet, although the altitude 
and atmospheric conditions are not described. This work also 
does not note any distinctive differences between the T-33 con- 
trails and the contrails of the other aircraft examined (F-86, 
F-100, B-47 , B-52) . 

• The 12,500-foot altitude is sufficiently different than the 1,000- 

foot altitude of Chevalier's tests that the local environment is 
not dominated by effects of ground heat transfer, airborne par- 
ticles, inversion layers, ground-generated pollutants, fog, and 
the like. 

Consequently, it is felt that early dissipation of the sort noted by 
Chevalier is not likely to occur with the T-33 tests. 

One interesting possibility is suggested by Chevalier's tests. If 
there is an unknown dissipative mechanism present as Chevalier's tests 
suggest, it would be most important to define that mechanism for possible 
consideration as an aircraft-borne dissipator device. ARC might) consider 
the possibility of performing lower altitude tests both in clear weather 
and in high humidity conditions (providing acceptable visibility and safety 
conditions) over Moffett Field where some atmospheric measurements can be 
made either by radiosonde or the Learjet instruments to attempt to dupli- 
cate Chevalier's results. 

/ 

VORTEX VELOCITY DATA REDUCTION 

‘ The Learjet will pass through the T-33 wake vortex system and obtain 

a record of the voltage outputs of the three wires of the probe along some 
path traversed by the probe. This section is concerned with the manner in 
which this data is analyzed to obtain information on the vortex characte- 
ristics. It is assumed that a horizontal traverse at right angles to the 
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wake will be made. Thus, the path of the probe will be a straight <pr 
curved line lying in a plane essentially normal to the vortex filament 
(denoted the crossflow plane). Further, it is assumed that some prior 
operation on the data will be done to convert the hot-wire voltages to 
velocities, to transform the three velocity components into a system having 
axes parallel to the principal axes of the probe aircraft, and to remove 
the velocity components at the hot wire due to aircraft motion. The re- 
sulting data will be velocity components normal to and coaxial with the 
path line followed by the aircraft in the crossflow plane and a third 
component normal to the crossflow plane (along the vortex axis) . Finally, 
it is assumed that time on the magnetic tape record and the aircraft velo- 
city components can be related such that two coordinates of the path of 
the hot-wire probe through the vortex system in an arbitrary coordinate 
system in the crossflow plane can be established. Thus, the velocity 
components at any point along this path line can be obtained from the data. 

The first part of this section contains a discussion of the types of 
vortex penetrations that can be achieved and the reconstruction of the 
vortex velocity system from the measured velocity components along the 
probe path. The second part then deals with the description and use of a 
computer program associated with the velocity reconstruction procedure. 

Only the two velocity components in the crossflow plane are considered in 
the following discussion. 

Methods for Reconstruction of the Vortex Characteristics 

Penetration of a vortex in a crossflow plane can be categorized as 
follows. When the probe flight path passes closer to the vortex center 
than its core radius, (r ; ) , a core penetration is said to result. A pass 
that misses the vortex center by more than the core radius is called a 
noncore penetration. In either case, the component of the vortex-induced 
perturbation velocity normal to the probe flight path, V , will exhibit a 
well-defined maximum and minimum. If a wing-tip vortex has circular, con- 
centric streamlines with negligible radial velocity, and if there are small 
influences from the opposite tip vortex, the peaks should be equal and 
opposite as shown in the sketch below. 
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Downwash induced by the opposite vortex would tend to shift the profile in 
the vertical direction so that asymmetry results. It should be noted that 
in the absence of any motion of the vortex center in the crossflow plane 
during the probe time, the component of the perturbation velocity along 
the probe flight path, V , should not change sign during the penetration. 

The resultant perturbation velocity is obtained from a vector addition 
of the components along and normal to the probe flight path. By erecting 
normals to the resultant velocity vectors, the vortex center can be 
located as shown in the sketch below. 



Probe 

flight 

path 
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From a knowledge of the vortex center and the resultant velocities, a curve 
of the variation of tangential velocity with radius can be determined. 

In fact, the tangential velocities V t are the resultant velocities V . 

The two types of vortex penetration can now be distinguished. For 
the case of a core penetration, the results will assume the appearance of 
the curve shown below. In this sketch the tangential velocity is plotted 
versus radius. Positive values of have an upward component, and 

negative values a downward component. No data points occur for values of 
the radius less than that for closest approach to the center. One sees 
well-defined maximum and minimum peaks in the tangential velocity variation 
The radial location where these occur is termed r . The radius of closest 
approach must be less than r for a core penetration to occur and for a 
well-defined maximum to occur. 



The peaks in the vortex tangential velocity V^_ must be clearly indicated 
before the core size can be estimated. In other words, some data points 
must exist on either side of the maximum and the minimum. 

If the probe does not enter the core, the same form of resultant 
velocity and vortex center determination may be used to construct a curve 
of tangential velocity versus radius. However, this curve will not show 
the well-defined maximum and minimum velocities, as shown below, since 
data will not be available at or within the core radius. 
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The method for performing the analysis just described is as follows. 
The axis system in the crossflow plane is shown in figure 10, together 
with the positive sense of the normal and coaxial velocity components, V n 
and V . The origin and the direction of the positive sense of the Z,Y 
coordinate system are somewhat arbitrary, in accordance with the following 
considerations. Figure 10 shows a left-to-right traverse, looking along 
the vortex towards the generating aircraft. The coordinate system shows 
Z generally down and Y to the right. It is suggested that an origin 
of this coordinate system be fixed for each set of data from a pass through 
the wake. The location of the orig;in is completely arbitrary and serves 
only to permit assignment of two length coodinates to the data taken at 
every instant of time during passage through the wake by means of the air- 
craft accelerometer data. The Z and Y axes do not have to be truly 
vertical and horizontal but only approximately so. However, the general 
arrangement of Z down and Y to the right should be maintained for all 
passes, whether right- to-left or left-to-right, in order to provide con- 
sistent signs and senses of rotation of the vortices in the data reduction. 

Pertinent angles necessary for the determination of the slopes of the 
normals to the resultant velocity vectors V r at field points (Z,Y) are 
also shown in figure 10. Two such lines emanating from points (Z^,Y^) and 
(Zj,Yj) are seen to intersect at point (Z^j,Y^j). The expressions for the 
intersection point coordinates are given by 
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The slope of line i is given by and the slope of line j is given 

by ni£j. These slopes are determined from the following expressions: 



and (AY/ AZ) ^ is the slope of the flight path at the i'th point, obtained 
from the coordinates of points along the path adjacent to the i'th point. 


If there are N data points given along a probe flight path, the 
number of unique intersection points is equal to 1/2 N(N-l). Ideally, all 
intersection points should coincide. In actuality, a collection of points 
is generated and after exclusion of those points that are obviously. in 
error , a suitable averaging scheme can be called upon to generate the 
apparent vortex center from the set of coordinate pairs. In 

the ideal case, lines drawn from the apparent vortex center to the field 
points where the resultant velocity vectors are constructed will be normal 
to those vectors. In practice, these lines will be near normal to the 
velocity vectors and their lengths are approximately equal to the radial 
distances from the vortex center to the tangential velocity vectors. 


The magnitude of the total circulation of the vortex can be 

determined from a potential vortex model using the vortex tangential 
velocities and the radial distances from the apparent vortex center . The 
approach followed here makes use of the following assumptions: 


(1) The vortex can be characterized by circular, concentric stream- 
lines in the crossflow plane. 

The radial velocity in the crossflow plane is negligible. 

The downwash effect from the opposite vortex is small. 

The vortex center does not move appreciably during the duration 
of the traverse. 


(2) 

(3) 

(4) 



These assumptions appear reasonable based on the following considerations. 
The small "intersection areas" obtained from the vortex center analysis of 
the wind-tunnel hot-wire data indicate that any radial velocities that do 
exist are small compared to the tangential velocities. Visual observations 
of the vortex size indicated no change in core diameter within 100 chords 
of the wing, which also indicates small radial velocities. Since the velo- 
city induced by one vortex of the pair at the other is only 5 percent of 
the maximum tangential velocity of that other vortex, the noncircularity of 
this effect should be small in determining vortex centers and circulations. 
This effect can be included later if it is found to be significant. Finally, 
the time to traverse the vortex is less than 0.01 second so that the vortex 
displacement in this time increment should be negligible . 

Since the potential vortex velocity distribution joins the inner 
velocity model at a radius (r ) about 2.5 times the core radius (fig. 2) , 
the following expression can be used for the calculation of a set of 
circulation values for those points having radii greater than 2.5 r x 

T = 2irr . V (20) 

°i i 

The average over the set of r oi ’s yields a value for the total circula- 
tion or vortex strength r'. 

Trailing Vortex Crossflow Program 

The purpose of this section is to describe the program which reduces 
the trailing vortex velocity data taken in some crossflow plane to a tabu- 
lation of vortex tangential velocities as a function of radial distance 
from the apparent vortex center. The program also computes an average 
total circulation associated with the vortex in the crossflow plane. 

Fundamentally, the program is based on the assumption that the vortex 
possesses circular, concentric streamlines with negligible radial velocity. 
The effects of the opposite tip vortex and vortex center motion are not 
accounted for. 

In the following, the program and its input and output will be 
described in detail. At the end of this section, a sample case illus- 
trates how the program handles a segment of data recently taken in the 
40- by 80-Foot Wind Tunnel at the Ames Research Center. 


42 



Calculation procedure .- The computer program proceeds through various 
stages as follows. After run identification and input of certain variables 
a set of field point coordinates and vortex perturbation velocity compo- 
nents coaxial and normal to the probe flight path in the crossflow plane 
are read in for one vortex traverse. The maximum and minimum values of 
the normal component and the corresponding field point coordinates are 
determined from the set of supplied data for later use. The program prints 
all of the input data and prints the peak values of the normal component 
velocity encountered during the traverse under consideration. 

The program proceeds to compute resultant velocities (V ) at all the 
input data points using the coaxial and normal velocity components. Since 
the induced velocities become very small at large distances from the vortex 
and can cause large errors in calculating the apparent vortex center, 
velocities smaller than a "cutoff" velocity are dropped. The cutoff 
velocity is calculated from input quantities. The remaining data set 
consists of the first point along the path where the velocity becomes 
greater than the cutoff velocity, all data points thereafter along the 
probe path through the vortex until that point where the velocity has 
dropped off again to the cutoff value. If any of the data points within 
this band have identically zero velocity components (for example, if the 
data point is judged bad and the velocity components set equal to zero) , 
the program will treat the velocity components as zero and continue to 
process that point. . 

After printing the field point coordinates within the range defined 
by the first and last velocities greater than the cutoff, the program 
determines the set of coordinates of the points of intersection of lines 
normal to the resultant velocity vectors. Resultant velocity vectors with 
magnitudes identically equal to zero are excluded from this calculation. 
Intersection points with Y coordinates outside of the range spanned by 

the field points where the previously calculated peaks in normal velocity 

‘ •) 

component occur are also excluded from further processing. The apparent 
vortex center is then computed from the subset of remaining intersection 
points by means of an averaging procedure. 

Upon completion of the determination of the coordinates of the 
apparent vortex center, the program continues with the calculation of the 
distances from the center to the location of the field points. Next, 
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these radial distances are tabulated in the output together with the 
corresponding resultant velocities. These velocities are in fact the 
vortex tangential velocities. The program then searches for the maximum 
and minimum velocity values and the corresponding radial distances. A 
check for core or noncore penetration follows. This is done by checking 
on the sense of the direction of the vortex tangential velocity at the 
field point immediately adjacent to the maximum or minimum peak towards 
the vortex center. 

Lastly, the program computes an average value for the total circu- 
lation from the set of circulation values obtained from equation (20) 
applied to all the tabulated vortex tangential velocities and corresponding 
radii. The program is arranged to handle cases consecutively . 

Program operation .- The trailing vortex crossflow program is written 
in Fortran IV language and has been run on the IBM 360/67 TSS batch system 
at Ames Research Center. With this machine, 500 data points can be pro- 
cessed although this number is not necessarily a maximum. If the program 
is to be run on a different computer with smaller storage, dimension state- 
ments need to be changed to permit operation on that machine . 

No tapes are required other than the standard input and output. The 
running time on the TSS batch system required to process about 300 data 
points is approximately 1.5 minutes. 

Program limitations and precautions .- Some care must be taken in 
specifying the variable FACTOR in the input. This positive variable multi- 
plied with the specified nominal value of the axial velocity UAX serves in 
the program as a lower bound, VMIN, for the magnitudes of the first and 
last resultant velocities in the crossflow planes treated by the program. 
When this lower bound is too low, the number of admissible vortex tangen- 
tial velocities becomes unduly large in its relation to the points of 
intersection determined by the program. The program is dimen- 

sioned to handle a set of 100 vortex tangential velocities with the first 
and last magnitudes larger than VMIN. If more than 100 are treated, the 
program will be stopped and the following message appears in the output: 
"DIMENSION OF INTERSECTION POINT COORDINATES NOT SUFFICIENT". An increase 
in the input variables FACTOR or UAX will raise the lower bound (cutoff) 
value, VMIN, and alleviate the problem. Alternatively, the array size of 
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the variables ZC(Z^) and YC(Y^j) . in the dimension statement can be 
enlarged if the storage capacity of the computer is large enough. 

The program writes out the coordinates of all intersection points as 
an option. This output can be time and paper consuming if the number of 
admitted vortex tangential velocities is large. It is recommended that 
the program be run with a given set of data first without the optional 
output. If the value for the total number of computed intersection points 
specified in the output is less than 200, the program can be rerun with the 
inclusion of the optional output if so desired. 

Description of input .- This section describes the input for the 
trailing vortex crossflow program. In the following discussion, the 
content of all input cards is described and, where appropriate, instruc- 
tions on generating the input quantities are given. All input variables 
are summarized at the end of this section in the order of appearance in 
the input deck. The input format for all cards is shown in figure 11; 
the item numbers refer to figure 11. 

The TSS 360 computer system at Ames Research Center offers an alter- 
nate way for reading in the set of field point coordinates with the two 
vortex perturbation velocity components. These data can be generated and 
stored in a data set. The trailing vortex crossflow program can be modi- 
fied to read from a data set with a minimum of alterations. This method 
is advantageous when large amounts of data are generated. 

Item 1:. The first card serves as identification and may contain any 
alphameric information required. This information is printed on the first 
page of the output. 

Item 2; The second card contains an axial velocity UAX and the axial 
location of the crossflow plane under consideration. The axial velocity 
can be tunnel velocity for wind-tunnel data or the speed of the vortex 
generating aircraft for flight data. This quantity is used only to cal- 
culate the cutoff velocity VMIN. In addition, this card specifies the 
fraction FACTOR of the axial velocity used to calculate the cutoff velo- 
city, the total number of data points supplied (500 maximum) and a vari- 
able controlling the optional output of the intersection point coordinates. 
As was noted above in the discussion concerning program limitations and 
precautions, the fraction specification multiplied by the nominal axial 
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velocity generates the lower bound or cutoff value, VMIN, for the magni- 
tudes of the first and last resultant velocities to be processed by the 
program. 

The choice of the quantity FACTOR (and therefore VMIN) is motivated 
by the desired number of data points to be included in the analysis and 
perhaps by the signal-to-noise ratio. As a rule, VMIN must be high enough 
so that the number of admitted data points does not exceed 100. At the 
same time, VMIN must be low enough to allow enough data points for an 
accurate construction of the vortex-induced flow field in the crossflow 
plane. In addition, if the signal-to-noise ratio is high, one might want 
to use a value for VMIN a certain amount above the velocity corresponding 
to the noise voltage level to improve the accuracy of the calculation. For 
the 40- by 80-Foot Wind Tunnel tests, FACTOR was found to lie in the range 
of 0.20 to 0.30. The tunnel velocity UAX was 90 feet per second so that 
VMIN ranged from 18 to 27 feet per second. Variable IZCYC controls 
optional output of intersection point coordinates. 

Item 3: This set of cards contains the field point coordinates and 

the perturbation velocity components induced at those points by the vortex 
in the crossflow plane. Each card contains a sequence number, the Z and 
Y coordinates and the velocity components coaxial and normal to the probe 
flight path. The sequential order of this input data must be the same as 
the order in which the experimental data was obtained. The coordinate 
system and the positive sense of the velocity components are shown in 
figure 10. It should be noted that any contribution to these velocity 
components by mechanisms other than the vortex should be excluded. The 
total number of cards is NFLDP (NFLDP <_ 500, Item 2) . 

A sample input deck is shown in figure 12. This data is part of a 
traverse performed through a vortex in the 40- by 80-Foot Wind Tunnel and 
serves as input to the sample case described below. In this input, the 
variable IZCYC is specified to be 1; consequently, the output will include 
the coordinates of all intersection points. 
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Program. 

Variable 


Comment 


Algebraic Symbol 
(if applicable) 

Item 1 


Any alphameric information may 
be put on this card for identi- 
fication of the calculation. 


Item 2 
UAX 

XSTAT 

FACTOR 

NFLDP 

IZCYC 

Item 3 

IN, 1,2,. . .NFLDP, 

ZFP(IF), Z (1) ,Z (2) , . . .Z (NFLDP) , 

YFP(IF), Y(l) ,Y(2) ,...Y (NFLDP) , 

VTANFP (IF) , V s (1) ,V g (2) , . . .V s (NFLDP) , 

VNORFP ( IF) V n (l) ,V n (2) ,. ..V n (NFLDP) 


Nominal axial velocity (tunnel 
speed) ; positive; feet per 
second 

Axial location of crossflow plane 
where data is obtained; feet 

Arbitrary positive fraction of 
UAX. NOTE : FACTOR x UAX = VMIN 

Total number of supplied data 
points; 1 < NFLDP <_. 500 

l 

IZCYC =0: coordinates of inter- 

section points not in output 
IZCYC = 1: coordinates of inter 

section points in output 


Data point coordinates, velocity 
components along and normal to 
the probe flight path, see fi- 
gure 8 for coordinate system; 
feet, feet per second. 

1 < IF <_ NFLDP. This input 
requires NFLDP cards. 


Description of output .- This section describes the output from the 
trailing vortex crossflow program. In the following, all items of output 
are enumerated and the contents briefly described. Sample output is shown 
in figure 13. 

The first pages identify the run and contain all the input data in 
connection with the probe traverse in the crossflow plane. The peaks in 
the normal component velocity and the associated field point coordinates 
determined by the program are also specified in this part of the output. 

In the ideal case, the peak values should be equal and opposite. 
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The next item, printed on a new page, specifies the lower bound or 
cutoff value VMIN for the magnitudes of the first and last resultant 
velocity vectors of the set to be treated further by the program. This 
is followed by a tabulation of the resultant velocities and coordinates 
starting with the first velocity magnitude larger than the cutoff value, 
then proceeding along the path through the vortex until the last velocity 
magnitude larger than VMIN. The first column is the index IF associated 
with the supplied data; the second column IW is another running index which 
starts with the value one. Index IW enumerates the set of field points' 
with the first and the last having resultant velocity magnitudes larger 
than VMIN. Some of the resultant velocity magnitudes within this band may 
be less than VMIN. The last three columns represent angles necessary for 

the geometric combination of lines normal to the resultant velocity vector. 

✓ 

In addition, the first and last field points of the set of admitted resul- 
tant velocities or vortex tangential velocities are also printed on this 
page. 

Optional output of coordinates appears on the next page. These coordi- 
nates locate the points of intersection of the lines erected normal to the 
resultant velocity vectors at the field points. This output can be lengthy 
and is therefore made optional through the control variable IZCYC in Item 2 
of the input. The first two columns are indices of the supplied data 
points, the second two columns are indices associated with the admissible 
resultant velocity vectors. The last two columns indicate the coordinates 
of the intersection of the two normals identified by the indices. 

On the next page, the total number of intersection points determined 
by the program is listed together with the number of admissible inter- 
section points. These are points with Y coordinates falling in the range 
of Y spanned by the peaks of the normal components of the vortex-induced 
velocities as identified at the end of the first pages of output. The 
coordinates of the apparent vortex center are also printed as the average 
values of the Z and Y coordinates of the subset of admitted inter- 
section points. 

The item of the output which appears on the next page is a tabulation 
of vortex tangential velocities (resultant velocities) and the distance 
from the apparent vortex center to the associated field points. This page 
also lists a statement concerning a check on core or noncore penetration 
of the traverse under consideration. 
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The maximum and minimum values of the vortex tangential velocities 
and the associated radii are printed on the next page. This information 
is printed for both core and noncore type penetrations. In either case, 
the peak values should be equal and opposite and the two radii equal for 
the ideal case of a single vortex with concentric, circular streamlines. 

If the traverse resulted in a core penetration, the average of the peak 
magnitudes equals the maximum vortex tangential velocity v x and the 
average of the two radii is r r . A tabulation follows of vortex strengths 
computed according to potential flow theory from the vortex tangential 
velocities and the radii. Only those velocities at radii larger than 
2.50 times the radius (r L ) where the peaks in velocity occur are included 
in the vortex strength computation. The last item of output is the average 
vortex strength or circulation determined from the tabulated values. 

Program listing .- The trailing vortex crossflow program is written 
in Fortran IV language for the IBM 360/67 computer. The program consists 
of a main program without subroutines. 
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PROGRAM TVCRS( INPUT .OUTPUT, TAPE 5 = 1 NPUT »TAPE6=0UTPUT ) 

PROGRAM TVCRS 

THIS PROGRAM COMPUTES THE STRUCTURE OF A TURBULENT VORTEX IN THE 
CROSSFLOW PLANE FROM EXPERIMENTAL FLOW FIELD DATA IN THE CROSSFLOW 
PLANE 

Y-AXIS TO THE RIGHT WHEN FACING FORWARD. Z-AXIS DOWN 

DIMENSION YFP< 500) ,ZFP( 500 ) ,VT ANFP ( 500 ) , VNQRFP I 500 ) . TALK! 20 ) . 
1VR5LT(500) , SLOP ELI 500) . VVRT AN( 500) 

DIMENSION ZC (100,100), YC 1100 ,100 ),RAD(10D) 

C 

C 

1 FORMAT (20A4) 

2 FORMAT (1H1.20A4//) 

3 FORMAT (3F10.5.2I5) 

4 FORMAT (///IX, 14, IX »57HF IELD POINTS ARE SPECIFIED IN THE CROSS FLO 
1W PLANE AT X= ,F10.5,1X,6H FT. / AX .25HN0M I NAL AXIAL VELOCITY IS, 
24X,F10.5,1X,6HFT/SEC///) 

5 FORMAT (///IX, 88HTH1S ANALYSIS IS CONCERNED WITH RESULTANT VELOCI 
1TIE5 IN THE CROSS FLOW PLANE LARGER THAN ,3X , FI 0 . 5, 1 X , 9HTIMES THE, 
2/1X , 14HAXI AL VELOCITY///) 

6 FORMAT ( I5»4F10 .5) 

7 FORMAT (///IX, 84HIN THIS VORTEX PENETRATION, THE PERTURBATION VELD 
1CITY NORMAL TO THE PROBE FLIGHT PATH/lX , 13HHAS A MAX IMUM.5X , F10 .5, 
21X,6HFT/SEC,2X,2HAT,2X,2HZ=,5X,F10.5,1X,6H FT. /42X ,2HY = , 5X ,F 10 .5 
3.1X.6H FT. /IX , 13H AND A MI N 1MUM ,5X , F10 .5, 1 X ,6HFT/ SEC , 2X ,2HAT, 2X , 

3 2HZ= 

4,5X,F10.5,1X,6H FT. /42X,2HY= , 5X, F10 . 5, IX , 6H FT. ///) 

8 FORMAT (///1X,7X,20HCOORDINATES IN CROSS , 1 3X .34HCR0SS FLOW PERTURB 
1 AT 1 ON VEL0CITIE5/12X.13HFL0W PLANE , FT, 19X.22HFR0M EXPER IMENT , FT/SE 
2C//4X.2HIF.3X, 

2 1HZ,14X,1HY,17X»7H CO AXIAL, 13X,6H NORMAL/) 

9 FORMAT (I5,F10.5,5X,F10.5,10X,F10.5,10X,F10.5) 

10 FORMAT ( 1H1 , 62HRESUL TAN T VELOCITIES WITH THE FIRST AND LAST LARGER 
1 THAN VM IN = , 

1 2X.F10.5, 

11X, 51HFT/S EC AT THE INDICATED FIELD POINT COORDINATES AND/ IX, 

1 46HSL 

2 OPES OF LINES NORMAL TO THE VELOCITY VEC TORS //3X ,2H I F , BX ,2H I W ,8X , 
31HZ , 15X , 1HY , AX, 

3 19HRESUL T . VEL., FT/SEC, 

3 2X.11HTHETAS, DEG., AX, llHTHETAL.DE 
AG. , AX.6HSL0PEL/ ) 

11 FORMAT (IX, IA,6X,IA,5X,F10.5,5X,F10.5,5X,F10.5,5X,F10.5,5X,F10.5, 
15X.F10.5) 

12 FORMAT (//1X.50HTHIS FLOWFIELD ACCOUNT STARTS AT DATA POINT NUMBER 
1,6X,3HIF=, I5/1X.52HWITH FIELD POINT COORDINATES IN THE CROSS FLOW 
2PLANE,5X,2HZ=,F10.5,2X,6H FT. /58X ,2HY= ,F10 .5.2X.6H FT. /1X.5AHT 
3 HE LAST OATA POINT NUMBER INCLUDED IN THIS ACCOUNT I S, 2X ,3H I F= , 1 5, 
A/1 X, 52HWITH FIELD POINT COORDINATES IN THE CROSS FLOW PLANE, 5X, 
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52HI * ,F10.5,2X,6H FT. /58X , 2HY* , F10 .5,2X ,6H FT. //) 

13 FORMAT (1H1, 

1 3X tlOOHCOOROl NATES OF THE POINTS OF INTERSECTION OF THE 

1LINES CONSTRUCTED AT POINTS (ZFP.YFP) NORMAL TO THE/ 

1 10X,42H RESULTANT VELO 

2CITY VECTORS AT THOSE PO INTS//9X ,2H I L, 8X ,2H JL, 8X ,2H I W, 8X ,2H JW, 13X , 
36HZC»FT .,9X,6HYC,FT./) 

14 FORMAT (9X,I3,7X,13,7X,I3,7X,13,8X,E12.5,5X,E12.5> 

15 FORMAT (1H1*4X»44HT0TAL NUMBER OF COMPUTED INTERSECTION POINTS, 26X, 

1 I5/5X, 40 HN UMBER OF ADMISSABLE INTERSECTION POINTS, 30X, I 5/5X , 54HAVE 
2RAGE Z COORDINATE OF ADMISSABLE INTERSECTION POINTS ,6X ,F10.5,1X ,3H 
3 FT . /5X , 54HA VERAGE Y COORDINATE OF ADMISSABLE INTERSECTION POINTS, 

4 6X , F10 .5, IX ,3HFT .///) 

16 FORMAT (1H1 ,5X, 3 8H VORTEX TANGENTIAL VELOCITIES AND RADI I//9X,2H1V, 

1 10X , 10HRAD IUS, FT .,5X,27HTANGENTIAL VELOC ITY , FT/SEC ./ ) 

17 FORMAT (8X,I3,9X,F10.5,15X,F10.5) 

18 FORMAT { ///5X,33HV0RTEX STR ENGTH ,GAMMAO , FT**2/$EC.//9X , 2HI V , 10X , 

1 10HRADI US, FT., 5X .27HTANGENT I AL VELOC ITY , FT/SEC ., 7X , 16HGAMMA0 , FT**2 
2 /SEC/) 

19 FORMAT (8X, I 3,9X ,F10 .5 , 1 5X , F 10 . 5 ,15X,F10.5 ) 

20 FORMAT <1H1,53HIN THIS VORTEX ACCOUNT, THE VORTEX TANGENTIAL VELOCI 
ITY /1X.13HHA5 ONE PEAK , 5X,F10*5, IX, 6H FT/SEC, 2X ,2HAT , 2X ,4HRA0 S , Fl 
20.5»2X»3HFT ./IX ,13 HAND ANOTHER , 5X , F10 .5, IX ,6HFT/SEC,2X ,2HAT, 2X ,4 
2HRAD=,F10.5,2X,3HFT .///) 

21 FORMAT (///5X.35HA VERAGE CIRCULATION IN OUTER REGION ,5X , FlO . 5, IX , 
19H FT** 2/SEC//) 

22 FORMAT I /// 1X.44HTH1S TRAVERSE RESULTED IN A CORE PENETRATION///) 

23 FORMAT ( ///l X, 48HTH I S TRAVERSE RESULTED IN A NON-CORE PENETRATION/ 
1//) 

24 FORMAT (1H1 » 5X , 59HD I MENS ION OF INTERSECTION POINTS COORDINATES NOT 
1 SUFFICIENT///) 

25 FORMAT ( /// IX , 31HAB0VE TABULATION HAS STOPPED AT ,2X , 3H I F = , 1 5 /// ) 

C 

1000 READ (5,1) TALK 
HRITE (6,2) TALK 
C 

C SPECIFY AXIAL FLOW VELOCITY UAX, AXIAL LOCATION, XSTAT , NUMBER OF 

C FIELD POINTS NFLDP . 

C FACTOR IS AN ARBITRARY FRACTION OF UAX AND 

C IS DETERMINED BY THE SIGNAL NOISE LEVEL 

C IZC YC s O , INTERSECTION POINTS NOT IN OUTPUT 

C IZC YC= 1 , INTERSECTION POINTS IN OUTPUT 

C 

READ (5,3) UAX, XSTAT, FACTOR, NFLDP, IZCYC 
VMIN=(FACTOR)*UAX 

C : •• • 

C SPECIFY COORDINATES IZFP,YFP) OF THE FIELD POINTS AND THE PERTUR- 

C BAT I ON VELOCITIES TANGENTI AL AND NORMAL TO THE PROBE FL IGHT PATH 

C IN THE CROSS FLOH PLANE, VTANFP AND VNORFP, FT AND FT/SEC 

C VELOCITIES POSITIVE HHEN COMPONENTS ALONG Z AND Y AXIS ARE POSITIVE 

C 

DO 102 1F*1, NFLDP 

102 READ (5,6) IN, ZFP( I F ) , YFP( 1 F ) , VTANFP ( I F ) , VNORFP ( IF ) 
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WRITE (6,4) NFLDP.XSTAT ,UAX 
WRITE (6,5) FACTOR 

VNQRMX=-1.0E-10 

VNORHN* 1.0E+10 

DO 131 I F=1 , NFLDP 

VNDRMX = AM AX 1 ( VNORMX , VNDRFP ( IF) ) 

131 VNQRMN- AMINK VNORHN, VNORFP ( IF) ) 

WRITE (6,8) 

DO 103 IF-1, NFLDP 

103 WRITE (6,9) IF , ZFP ( I F) , YFP ( I F) , VTANFP{ I F ) , VNQR FP ( I F) 

LOCATE COORDINATES OF DATA POINTS WHERE MAXIMUM AND MINIMUM NORMAL 
PERTURBATION VELOCITIES OCCUR 

DO 107 I F=1 , NFLDP 
DIFFMX*VNORMX~VNQRFP ( IF) 

DI FFMN* VNORFP ( I F )-VNORMN 

IF (OIFFMX.LT.1.0E-07) IFMAX*IF 

IF (DIFFMN.LT. 1.0E-07) IFMIN=IF < 

107 CONTINUE 

WRITE (6,7) VNORMX, ZFP ( I FMAX ) , Y FP( I FMAX ) , VNORHN, ZFP(IFNIN), YFP 
KIFMIN) 

DIFFY*YFP(IFMAX)-YFP(IFMIN) 


COMPUTE RESULTANT VELOCITIES AT POINTS (ZFP, YFP) IN THE CROSSFLOW 
PLANE AT X*XST AT ION • IF MAGNITUDE IS LARGER THAN VH1N, CALCULATE 
SLOPE OF THE NORMAL TO THE RESULTANT VELOCITY VECTOR ANO PRINT 
THEM AS A SUBSET OF THE SUPPLIED DATA. 

WRITE (6,10) VMIN 
IF I RST = 0 

iw = 0 

DO 100 IF=1, NFLDP 

VRSLT(IF)*SQRT(VTANFP< IF)**2*VN0RFP( IF)**2) 

IF (IFIRST.EQ.O. AND. VRSL T ( I F ) .GT . VMIN) GO TO 104 
IF ( IFIRST.GT.O.AND.VRSLT( IFJ.LT.1.0E-07) GO TO 941 
IF ( IFIRST.GT.O.AND. 

1 1FMAX.GT. IFNIN.AND.IF.GTaIFHIN.AND.IF.LT.! FMAX . AND .VRSL T (IF). 
1LT.VMIN) GO TO 941 
IF (IFIRST.GT.O.AND. 

1 IFMI N.GT . I FMAX. AND. IF. GT .IFMAX.AND.1E.LT .IFMIN . AND .VRSL T (IF). 
1LT.VMIN) GO TO 941 

IF ( IF IRST.GT. O.AND . IFMI N.GT. I FMAX. AND . I F.GT .1 FM IN . AND • VRSL T ( I F ) . 
1LT.VM1N) GO TO 942 

IF (IFIRST.GT.O.AND . IFHAX.GT .IFMIN. AND. IF. GT .1 FM AX . AND .VRSLT ( I F ) . 
1LT.VM1N) GO TO 942 
IF ( 1FIRST.GT.0) GO TO 101 
GO TO 100 
104 IF IRST*I F 
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GO TO 101 

941 VRSLT( IF)=0.0 
101 IW= I W+l 

IF ( VNORFP ( IF) . GT.O .0) VVRTANI IW)=VR5LT( IF) 

IF (VNQRFP( IF) .LT.0.0) VVRT AN( I W )=-VRSLT ( I F ) 

IFP1*IF*1 

DELTAY*YFP( IFP1)-YFP( IF) 

DELTAZ=ZFP(IFP1)-ZFP(IF) 

THETAS = ATAN2 ( DELTAY,DELTAZ) 

IF (THETAS.LT. 0.0) THETAS=6.28318531*THETAS 
ALPHAV=ATAN ( VTANFP ( 1 F ) /VNORFP ( I F) ) 

THETAL =THETAS*ALPHAV 
SL0PEL ( I W) = TAN (THETAL) 

THETAD=THETAS*57. 2957795 
THETLD^THET AL*57 .2957795 

WRITE (6,11) IF,IW,ZFP(IF), YFP ( IF),VVRTAN(IW ), THETA 0 ,THETLD , SLOP EL 
1 ( 1 W ) 

HAST® IF 
100 CONTINUE 
GO TO 943 

942 WRITE (6,25) IF 

943 CONTINUE 


WRITE (6,12) IFIRST,ZFP( IFIRST) ,YFP( IFIRST) , HAST, ZFP( I LAST ) , 

1YFP ( 1LAST) 

DETERMINE COORDINATES OF THE POINTS OF INTERSECTION OF THE LINES 
CONSTRUCTED AT POINTS (ZFP.YFP) NORMAL TO THE RESULTANT VELOCITY 
VECTORS AT THOSE POINTS 

ILSTM1= ILAST-1 

DO 105 I L= IFIRST, ILSTM1 

JIS=IL*1 

DO 105 JL= JLS, 1 LAST 
1W-IL-IFIRST+1 
JW=JL— IFIRST+1 

IF ( IW.GT.IOO.OR.JW.GT.IOO) GO TO 196 
GO TO 197 

196 WRITE (6,24) 

STOP 

197 CONTINUE 

IF ( ABS ( VVRTANI IWJJ.LT.l .OE-07 .0R.A6S(VVRTAN(JW) ).LT .1.0E-Q7) GO 
1T0 124 
GO TO 106 

124 ZC( IW, JW)«1.0E-12 
YC( IW, JW)*1.0E-12 
GO TO 105 

106 IF ( ABS ( SLOP EL ( JW) -SLOPELl IW I ) .LT.1.0E-07) GO TO 127 

ZC( IW, JW)«-( (YFP(IL)-SLOPEL( 1W)*ZFP< IL))-( YFP( JL I-SLOPEL ( JW ) * 

1ZFP ( JL ) ) )/ 

i ( slope l ( iw ) - slop el ( jw) ) 



YC( IW,JW)= YFP( IL)-SLOPEL( IN )•( ZFP ( I L) -ZC( IW,JW) > 

GO TO 105 

127 ZC ( I H , JW ) = 1 .0E+12 
YCC IW,JW)=1 .OE+12 
105 CONTINUE 

C 

C PRINT COORDINATES OF ALL POINTS OF INTERSECTION t GO 00 OR BAD 
C 

IF ( IZCYC.EQ.OI GO TO 198 
WRITE (6,13) 

DO 126 IL = IFIRST , I L5TM 1 
JLS s IL + 1 

DO 126 JL- JLS, ILAST 
IW*IL-IF IRST+1 
JW« JL-IFIRST+1 

126 WRITE (6,19) I L , JL , I W, JW ,ZC < IW , JW) , YC( I W , JW ) 

198 CONTINUE 
C 

C EXCLUDE INTERSECTION POINTS WITH YC OUTSIDE THE REGION SPANNED BY 

C YFPUFMAX) AND YFP(IFNIN) 

C 

IP I B=0 

DO 108 IL*IFIRST,1LSTH1 
JLS = IL + 1 

DO 106 JL- JLS, 1 LAST 
IW- I L—I FIR5T+1 
JW*JL-IF1RST+1 
IF (DIFFY.LT. 0.0) GO TO 110 
GO TO 111 

110 IF (YC(IW,JW).LT. YFPUFMAX). OR .YC ( I W, JW ) .GT .YFP ( 1 FN1 N ) ) GO TO 113 
GO TO 108 

111 IF ( YC (IW.JW ).GT. YFPUFMAX) . OR .YC ( I W , JW ) .LT .YFP ( I FM IN) ) GO TO 119 
GO TO 108 

113 ZC( IW,JW)*0.0 
YC ( I W, JW )=0 .0 
IP IB *IP I B+l 
GO TO 115 
119 ZC( IW, JW1-0.0 
YC( IW, JW)=0.0 
IP I B = I P 18+1 
115 IP IBT» I P IB 
108 CONTINUE 
C 

C CALCULATE AVERAGE COORDINATES ( Z CAVE ,YC AVE ) OF GOOD INTERSECTION 

C POINTS 
C 

ZCA VE»0 .0 
YCAVE*0 .0 

JLT0T»(ILAST-IFIRST)+1 
NIPTOT= ( JLTOT*( JLTOT-1 ) ) /2 
NIPG*NlPTOT-IPIBT 
AlfIPG=NIPG 

DO 116 1L»1FIRST,1LSTM1 
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ULS = IL+1 

DO 116 JL* JLS, I LAST 
1W= I L- 1 F IRST+1 
JW= JL-I FIRST+1 
ZCAVE=ZCAVE+ZC ( IM,JW)/ANIPG 
YCAVE* YCAVE+YC ( IW,JW)/ANIPG 

116 CONTINUE 
C 

C PRINT COORDINATES ( ZCAVE, YCAVE ) OF THE AVERAGE INTERSECTION POINT 

C ALSO THE TOTAL NUMBER OF INTERSECTION POINTS AND THE NUMBER OF 

C GOOD INTERSECTION POINTS 

C 

WRITE (6,15) NIP TOT , NIP G,ZCAVE,Y CAVE 

C 

C COMPUTE RADII OF THE VORTEX AS THE DISTANCE FROM THE AVERAGE 

C INTERSECTION POINT (VORTEX CENTER) TO THE DATA POINTS 

C 

WRITE (6,16) 

DO 117 I F= I F IRST , I LAST 
I V* I F-I F IR ST-*1 

RAD ( I V) = SQRT((ZFP(1F)-ZCAVE)**2-MYFP(IF)-YCAVE)**2) 

IF( ABS (VVRTAN(lV)J .LT • 1 .0E-07) RAD(IV|*0.0 

117 WRITE (6,17) IV,RAO(IV),VVRTAN(IV) 

C 

c 

VTANMX*-1.0E-1Q 
VTANMN=1 .0E + 10 
DD 141 I F = I F IRST , ILAST 
IV* I F— 1 F IRS T+l 

VTANMX * AHAX 1 ( VT ANMX , VVRTAN ( IV)) 

141 VTANHN=AHIN1(VT ANMN , VVRTAN (IV)) 

DO 142 IF=IFIRST, ILAST 

IV* I F-IF 1R5T+1 
DIFFMX=VTANMX-VVRTAN(IV) 

DI FFMN = VVR T AN( 1 V)-VTANMN 
IF(DIFFMX.LT.l .0E-07 ) IVMAX=IV 
IF (DIFFMN.LT. 1.0E-07) IVMIN*IV 

142 CONTINUE 
C 

C CHECK FOR CORE PENETRATION 

C 

IF ( IVHAX.LT. IVMIN) GD TO 192 
IVMNP1 * I VMIN+1 

IF ( VVRTANt IVMNPl) .LT.O.O) WRITE (6,22) 

IF ( VVRTAN ( l VMNP1) .GT.O. 0) WRITE (6,23) 

GO TO 193 

192 IVHXP1*IVHAX+1 

IF ( VVR TAN ( I VHXP1) . GT.O .0) WRITE(6,22) 

IF ( VVRTANt IVMXP1) .LT.O.O) WRITE(6,23) 

193 CONTINUE 
C 

WRITE (6,20) VTANMX, RAD( IVMAX) ,VTANMN,RAD( IVMIN) 



C CALCULATE CIRCULATION IN OUTER REGION (POTENTIAL VORTEX) 
C 

RADlMRADt I VMAXl+RADI I VN IN) 1/2.0 

IVNG=0 

5UMGAM=0 .0 

WRITE (6,18) 

RADMN=2.50*RAD1 
DO 118 IF=IFIRST,ILAST 
IV® I F-I F IRS T+l 

1F(RAD (I V) .LT.RADMN ) GO TO 119 
GO TO 120 

119 IVNG=I VNG+ 1 
VORTAN =0 .0 
GAHM A0=0 .0 
GO TO 118 

120 V0RTAN = AB5 ( VVR TAN( I V ) ) 

GANMA0=6.2831*RAD( IV)»VORTAN 
SUMGAM*SUMGAH+GAMHAO 

121 WRITE (6,19) IV,RAD( IV) , VORTAN, GANNAO 
118 CONTINUE 

ATOT=JLTOT- I VNG 
GAMA V®SUMGAM/ATQT 
WRITE (6,21) GAMAV 
GO TO 1000 
END 
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Sample Case 

The sample case concerns a segment of data obtained during one tra- 
verse through a wing tip vortex in a test performed recently in the 40- by 

80-Foot Wind Tunnel at Ames Research Center (Contract NAS2-6719) . The 
axial location of the crossflow plane was 55 chords behind a 6-inch chord, 
32-inch span rectangular wing at 12° angle of attack. A three-axis, hot- 
wire probe measured the axial, coaxial, and normal velocity components. 

The probe flight path was an arc with 7.33 feet radius. For purposes of 
this analysis the probe flight is assumed straight during the length of 
vortex penetration. 

The input data for this sample case is shown in figure 12; the output 

is listed in figure 13. A plot of the resultant velocities and the appa- 

rent vortex center is shown in figure 14. Using the data of figure 14, 
the magnitudes of the vortex tangential velocities are plotted as a function 
of radius in figure 15. This plot is a graphical representation of the 
corresponding tabulation in the program output of figure 13. 

The input for this sample case comprises 40 data points, and the 
normal component of the perturbation velocity has a maximum and minimum 
almost egual in magnitude, an indication of a traverse through a circular 
vortex with concentric streamlines. The number of data points admitted by 
the lower bound or cutoff velocity of 18.0 feet per second is 35. This 
particular traverse resulted in a core penetration with nearly equal peaks 
in vortex tangential velocity. The total circulation strength r Q of 
this vortex calculated by the program is 21.9 feet 2 per second. With the 
assumption of elliptical loading on the rectangular wing model, the 
theoretical circulation strength for the flow conditions for this test 
(a = 12°) was estimated to be 20.2 feet 2 per second, which is in good 
agreement with the program computed value. 

CONCLUDING REMARKS 

Several final observations should be made with respect to the planning 
and conduct of the flight tests to investigate wake vortices. First, it 
should be emphasized that there are essentially no reliable, comprehensive 
data on flight measurements of wake vortices. The ARC tests represent an 
excellent opportunity to make a really fundamental contribution to the 
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knowledge of wake vortex behavior. Although a primary goal of the tests 
is to measure fully developed vortex properties to assess the wake- 
encounter problem, it should also be a primary goal to obtain sufficient 
measurements to shed some light on the dissipation problem. Since the 
dissipation mechanisms are so poorly understood, as much auxiliary data 
on such effects as temperature, humidity, and turbulence should be taken 
as are consistent with the time and funds available. These data will 
serve many purposes to many investigators, and the thoroughness with which 
the test parameters and conditions are defined will be very important to 
users of the data. As an aid in planning the data collection and reduction 
effort, a summary of the recommended measurements, with some commentary and 
a reference to the discussion in this report, is presented as figure 16. 

NEAR would like to emphasize the importance of good preparation in 
using hot-wire anemometry in a program of this sort, based on our experi- 
ence in the 40- by 80-Foot Wind Tunnel tests. We had considerable diffi- 
culties with the mechanical integrity of the probe, the recording equipment, 
the calibration, the digitization process, and the data reduction procedure. 
The more experience that can be gained before the tests with the various 
facets of using the instrument and reducing the data, the more efficiently 
the data collection and reduction procedure will be. 

The data reduction scheme was developed making certain assumptions 
regarding the symmetry of a vortex and the lack of radial velocity and 
ignoring effects of the other vortex of the pair and any vortex motion 
induced during the probe by the presence of the probe itself. This was 
felt to be a logical place to start. The success with which the program 
handled the wind-tunnel vortex data is an indication that the method should 
work well for the flight case. When the flight data is analyzed, however, 
considerable attention should be given to the behavior of the results to 
insure that the assumptions are not being violated. 

Finally, NEAR would like to emphasize again the importance of 
obtaining core penetrations, since considerably more information on the 
structure of the vortex is obtained than if only noncore penetrations are 
achieved. 


Nielsen Engineering & Research, Inc. 
Mountain View, California 
January, 1973 
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TABLE I 


PERTINENT CHARACTERISTICS OF T-33 AIRCRAFT 


Weight 

(at time vortices are laid down) 

11,750 pounds 

Wing Span 

37.54 feet 

Wing Taper Ratio 

2.63 

Wing Aspect Ratio 

6.0 

Wing Quarter Chord Sweep 

0 degree 

Wing Twist 

1-1/2 degrees (washout) 

Wing Airfoil Section 

NACA 65-213, a = 0.5 


TABLE II 

PROPERTIES OF THE NACA STANDARD ATMOSPHERE 
AT 12,500 FEET ALTITUDE 


Density 

Temperature 

Pressure 

Kinematic Viscosity 


p = 0.0.01621 slugs per feet 3 
T = 14.5° F 
p = 1320 psfa 

v = 2.14x10 4 feet 2 per second 


Speed of Sound 


c = 1067.4 feet per second 



TABLE III 


LIST. OF CONDITIONS USED IN 
CONDENSATION CALCULATIONS 



DCH-2 

T-33 

B-47E 

Gross Weight (lbs) 

4,100 

11,750 

200,000 

Span (ft) 

48 

37.5 

116 

Speed (ft/sec) 

150 

300 

880 

Altitude (ft) 

1,000 

12,500 

35,000 

r Q (ft 2 /sec) 

314 

819.5 

3,390 

r s (ft) 

1.1 

0.86 

2.66 

T l (ft 2 /sec) 

195 

509 

2 ,100 

v (ft/sec) 

28.2 

94 

126 



■TABLE IV 

RESULTS OF CONDENSATION CALCULATIONS 



DCH-2 

T-33 

B-47E 

Poo (ps f a)' 

2,040 

1,317 

498 

P M ~ P c (P sf ) 

5.8 

45.5 

37.1 

T (° F) 

OO 

55.4 

14.4 

-65.8 

fa 

o 

O 

Eh 

1 

8 

Eh 

0.3 

3.5 

6.3 

Pounds of water vapor 
condensed in core per 
pound of mixture* 

(w - 0) ) 

' oo c 

3x10" 5 

1.4xl0" 4 

1.3x10" 

Percentage of ambient 
water vapor condensed 




in core* 




/(i) - OJ \ 




C .. > io ° 

0.3% 

4.6% 

13% 

Temperature rise in core 
due to condensation* (° F) 

0.1 

0.6 

0.06 

Ratio of momentum imparted 
to condensed water vapor 
to momentum of core flow 

4.6xl0 -5 

2.1xl0" 4 

2xl0" 5 


*Based on the assumption that the ambient air is saturated at T^ 





(b) Wires 2 and 3 


Figure 1.- Example of a set of oscilloscope 
photographs of a core penetration with a 
hot-wire probe, taken in the 40- by 80 -Foot 
Wind Tunnel trailing vortex tests . 





turbulent vortex 
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by 10-Foot Wind Tunnel 





Figure 5.- Prediction of vortex deflection for flight conditions 



r 


Aircraft Stations: IZZtZ. Nose Cockpit Wing Root Chord 



x, feet 
(a) a = 0* 

Figure 6.- Perturbation velocity components around the forward 
portion of the Lear jet in the vertical plane of symmetry. 
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Figure 7.- Perturbation velocity components around th* 
forward portion of the Learjet in a vertical plane 
10.25 feet outboard of the vertical plane of 
symmetry (60 percent semispan) . 
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120 



(a) Induced angle of attack. 




Figure 8.- Angle of attack and vertical acceleration 
induced on Lear jet penetrating T-33 wake 
at 300 ft/sec. 




Figure 9.- Saturation pressure for water vapor 
at low temperatures. 
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Figure 10.- Vortex crossflow plane coordinate system. 















































SAMPLE CASE, CIRCULAR VORTEX, 40X80 TUNNEL DATA, RUN 16, FILE 38, RECORD 23 


90.0 

27.5 

O.i 

20 40 

1 

1 

0.0 

1.85625 

1.073 

15.862 

2 

0.0 

1.86750 

-0.205 

15.131 

.3 

0.0 

1.87875 

-1.170 

20.424 

4 

0.0 

1.890 

-0.368 

18.55899 

5 

0.0 

1.90125 

-1 .998 

15.901 

6 

0.0 

1.91250 

-0.498 

21.21999 

7 

0 .0 

1.92375 

-1 .880 

25.77399 

8 

0.0 

1.93500 

-3.678 

29.97499 

9 

0.0 

1.94625 

-3.105 

2 7.76799 

10 

0.0 

1.95750 

-3.830 

32.81299 

11 

3.0 

1.96875 

-5.043 

39.37999 

12 

0.0 

1.980 

-7.538 

40.894 

13 

0.0 

1.99125 

-11.493 

49.23799 

14 

0.0 

2.0025 

-17.69899 

45.72400 

15 

0.0 

2.01375 

-21 .508 

43.77599 

16 

0.0 

2.02500 

-28.700 

25.23099 

17 

0.0 

2.03625 

-32.310 

-2.357 

18 

0 .0 

2.0475 

-31 .06799 

-28.83299 

19 

0.0 

2.05875 

-23.642 

-43.63499 

20 

0.0 

2.070 

-17.53999 

-50.877 

21 

0.0 

2.08125 

-9.3640 

-51.750 

22 

0.0 

2.09250 

-8.116 

-47.42099 

23 

0.0 

2.10375 

-4.774 

-42.53799 

24 

0.0 

2.1150 

-4.549 

-36.670 

25 

0.0 

2.12625 

-0.417 

-34.06299 

26 

0.0 

2.13750 

-0.211 

-30.547 

27 

0.0 

2.14875 

-0.883 

-29.33199 

28 

0.0 

2.160 

-0.425 

-28.05399 

29 

0.0 

2.17125 

00.367 

-26.758 

30 

0.0 

2.1825 

1 .792 

-24.489 

31 

0.0 

2.19375 

-14.561 

-26.64099 

32 

0.0 

2.205 

4.087 

-22.920 

33 

0 .0 

2.21625 

2.483 

-21.39899 

34 

0.0 

2.2275 

0.866 

-19.51199 

35 

o.o 

2.23875 

2.807 

-18.92599 

36 

0.0 

2.25 

-0.234 

-18.713 

37 

0.0 

2.26125 

2 .2416 

-20.04599 

38 

0.0 

2.27250 

0.849 

-15.956 

39 

0.0 

2.28375 

-0.683 

-15.74 

40 

0.0 

2.29500 

-0 .965 

-13.279 


Figure 12.- Listing of sample input. 
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SAMPLE CASE, CIRCULAR VORTEX, 40X80 TUNNEL DATA, RUN 16, FILE 38, RECORD 23 


40 FIELD POINTS ARE SPECIFIED IN THE CROSS FLOH PLANE AT X« 27.50000 FT. 
NOMINAL AXIAL VELOCITY IS 90.00000 FT/SEC 


THIS ANALYSIS IS CONCERNED WITH RESULTANT VELOCITIES IN THE CROSS FLOW PLANE LARGER THAN .20000 TINES THE 

AXIAL VELOCITY 



COORDINATES IN 

CROSS 

CROSS FLOW 

PERTURBATION VELOCITIES 


FLOW PLANE 

.FT 

FROM EXPERIMENT, FT/SEC 

IF 

Z 

Y 

COAXIAL 

NORMAL 

1 

0.00000 

1.85625 

1 .07300 

15.86200 

2 

0.00000 

1 .86750 

-.20500 

15.13100 

3 

0.00000 

1.87875 

-1.17000 

20.42400 

4 

0.00000 

1.89000 

-.36800 

18.55899 

5 

0.00000 

1.90125 

-1 .99800 

15.90100 

6 

0.00000 

1 .91250 

-.49800 

21.21999 

7 

0.00000 

1.92375 

-1 .88000 

25.77399 

8 

0.00000 

1.93500 

-3.67800 

29.97499 

9 

0.00000 

1 .94625 

-3.10500 

27.76799 

10 

0.00000 

1.95750 

-3.83000 

32.81299 

11 

0.00000 

1.96875 

-5.04300 

39.37999 

12 

0.00000 

1 .98000 

-7.53800 

40.89400 

13 

0.00000 

1.99125 

-11.49300 

49.23 799 

14 

0.00000 

2.00250 

-17.69899 

45.72400 

15 

0.00000 

2.01375 

-21 .S0800 

43.77599 

16 

0.00000 

2.02500 

-28.70000 

25.23099 

17 

0.00000 

2.03625 

-32.31000 

-2.35700 

18 

0.00000 

2.04750 

-31.06799 

-28.83299 

19 

0.00000 

2.05875 

-23.64200 

-43.63499 

20 

0.00000 

2.07000 

-17.53999 

-50.87700 

21 

0.00000 

2.08125 

-9.36400 

-51.75000 

22 

0.00000 

2.09250 

-8.11600 

-47.42099 

23 

0.00000 

2.10375 

-4.77400 

-42.53799 

24 

0.00000 

2.11500 

-4.54900 

-36.67000 

25 

0.00000 

2.12625 

-.41700 

-34.06299 

26 

0.00090 

2.13750 

-.21100 

-30.54700 

27 

0.00000 

2.14875 

- .8 B 300 

-29.33199 

28 

0.00000 

2.16000 

-.42500 

-28.05399 

29 

0.00000 

2.17125 

.36703 

-26.75800 

30 

0.00000 

2.18250 

1.79200 

-24.48900 

31 

0.00000 

2.19375 

-14.56100 

-26.64099 

32 

0.00000 

2.20500 

4.08700 

-22.92000 

33 

0.00000 

2.21625 

2.48300 

-21 .39899 

34 

0.00000 

2.22750 

.86600 

-19.51199 

35 

0.00000 

2.23875 

2.80700 

-18.92599 

36 

0.00030 

2.25000 

-.23400 

-18.71300 

37 

0.00000 

2.26125 

2 .24160 

-20.04599 

38 

0.00000 

2.27250 

.84900 

-15.95600 

39 

0.00000 

2.28375 

-.68300 

-15.74000 

40 

0.00000 

2.29500 

-.96500 

-13.27900 


IN THIS VORTEX PENETR AT ION ,THE PERTURBATION VELOCITY NORMAL TO THE PROBE FLIGHT PATH 


HAS A 

MAXIMUM 

49.23799 

FT/SEC 

AT 

Z« 

Y» 

0.00000 

1.99125 

FT. 

FT. 

AND A 

MINIMUM 

-51.75000 

FT/SEC 

AT 

Z* 
Y « 

0.00000 

2.08125 

FT. 

FT. 


Figure 13.- Listing of sample output. 



RESULTANT VELOCITIES WITH THE FIRST AND LAST LARGER THAN VNIN = 18.00000 FT/SEC AT THE INDICATED FIELD POINT COORDINATES AND 

SLOPES OF LINES NORMAL TO THE VELOCITY VECTORS 
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RESULTANT VELOCITY VECTORS AT TWOSE POINTS 
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Figure 13.- Continued, 



789776-03 2 . 01683 E *00 9 26 7 29 1 . 29916 E -03 1 . 957386*00 


ooooooooo 0 ' 4 oooooe>oooooooooooooooo*-oooooooo©oooooooooooooooo- 90 ©o«>*ooooao 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooeooopooo 

UOJUJlUlUIUlilUUiUiliJUMUiuUjUJlUlUlUUJWUJUnWlUliJIillUUjUJUjIlJkUlliUliJUJUMUtUUiaiUUJlUtUIUIUlUUJiUUjUlujunuiUWlLlUUjlUUJlUWIUUJlULUUJ 

as r~ ** a> in j a) ® rg p* inr*j»r'cot* > o»9Hoj«nnm*j <t m m c o No* <oiftO(Mn^»ifidWMnrtwn 

o< « » >r h ^ r^r-*-«i>oooor»oooooooaooof>CT> otT'^^-o^r-^c^o ooooo ooooooooooooo.o.hxi^cd »»iaoooooo 




3 rg rg rg 
J O O O 
I I t 


I I I I ♦ » 

I IU l>J IU WBJ *U 
}m**O40u\ 

- o ® rg eg in 
) mO 1 » <3 

1 H <1 IM rt IB 

- mom-J 


m m 
O O 
I I 

OJ LU U 

Al .4 t 
ro O - 
O O' C 


m «n p 

0 o c 

1 I I 


J m 43 

r -a rg 

< o* *4 


i m <n 
> o o 
i i 
UJ UJ 
in — 

rg >4 

co <n 


Oi—Aji—r— — <r-rg 

3 m *0 r- «-• *4 r* 


rt rt nif) w i 

o o o o o < 
1(111 

|f / 111 111 iu iu 

NniMQ*" 
o .4 in ® r- 
O* O ® f- «- 
OBO'ON 


*4 «o 
r- ffl 
a- cn 


ONIA 
(M*lO 
*in<M 
I- ~ rg 
m r- m 


INN NO » 
3 O O O O < 

III* 

j IU UJ LU UJ I. 
P' O MA 

ro r*> «n 
nN ip o> 

rtf' o® 

* (M W O 


in o> m 
m -9 4> t*3 
m — • m «vi 
013 

O' in m o 


® ® 
<n m 

O' m 


UJ UJ UJ 
hNN 
« IMN 

to >r r~ 

Ift 

•o r- -« 


IU UJ 
« O' 

a *7 
« <n 


i ip ip pi pi ip ip m (A ip n 

soooooooooo 
i i l i l i i i i i 

UJ UJ UJ UJ UJ UJ UJ UJ UJ UJ 

O lA ^ co <0 fp| 00 P'- .— Pg 

ONNMnmaftnQ 
N>T«^*JIA(MAhQ 
imia n jjr* hia o® > 
IP-T"tfOONh«nlfl 


« *NJ 

0 o 

1 I 

UJ UJ I 

vn in * 
o o* . 
-9 >4 f 
-■ 43 f 


M -■ 
O O 
I I 


I I 1 I I I 


00 

1 

ii i ii i |u UJ cu UJ UJ UJ UJ iu 
oo m ~~ 

» m 
r- r* 
in m 

o r- 


m rg rg 

00 O 

1 I I 


i i«P*in JJ ho »on 


>*o 
r o> o 

> f» ® 


O 4 
CD 4 
H IA 

O' O' 


SPIN 
ON O 
00 O O 


^NfpiANho 3 NHip-<«r'r.r'®mO' 0 ' 0 '( 3 fflO-® 0 'HHaNiP>rN®ir'H t HPiasij 44 r'oco 09 ' 0 , co 4 ‘a 30 , HH^NiPipAj^mHNiPN< 3 r*»>^«HH 


in -o r» asoo— •pgm*^*nO' 0 ^«njfpi>ymor^(no>o— •fMm'Tkn 4naa.OHNm'JiAOHNmtfK\4f'BooHNm>3iA4p.®froHNiP'tfiA hnpm 

NNNNNiPionmnn H^HHMHHHMHNNNNNNNNNNmmnninmHHHHHHHHHHNNNNNNNNNNpwnnwwrtHHHHHH 




40 00 (D co co co c 


r>(ooo«H<Nin4<A 

NNNmnpMApnn 


ss 


p'jm'»'n^r-ooO'OH*pgm^*n'Or*® 0 ' 0 «HN»n^-ing)p^nj(n^»n^r^cD<^Q 

HHH^^HHHIMftJNN'MNNNNNmiPmOlPlPfAnHHHHHMHHN 


HNn 4 lfi«h» 9 iOHNi*) 4 lA 4 MA 4 IA 4 h® 

NNNNNNNNNIPlPlPIPlPIPlPIPHHHHHH 


O' O' O' O' O' O' O O' O' O' o> O QOOOOOOOOOOOOOOOOOOOOOOOQO 


rg rg 


rO O' 
O' r- 
-* O 

rg pg 
o o 


44 WlArt. 
>4 o ® 4) i>- . 
OHIAON' 

O O O o o c 


► ♦ ♦ ♦ ♦ 

J UJ UJ UJ UJ I 

> r- m O' in 

JOBhO 

3NN4N 
34 IP4« 

3 0 0 0 0* 


4 ® ® n 4 
N 4 M 4 (P 
NO>44IA 
4 CO 4 4 a 
O' O' O* ® eg 


4 o o o o < 
3 O O O O < 
•♦♦♦♦• 
j UJ Uj uj uj i 
o o r- o c 

O' 43 O 43 U 
„« O' — 0 

O O' H IA V 

H in r* & u 


NNNNNNN-MNNNNH^HHrtMNNPlHN- 


3-9-940 
> gf *4 ® r- 
«4(MflH 

- O' ^ o *r 
' <• co *4 r*> 

9 rg -9 — <g 


oooooooooooooooooooooooo 
oooooooooooooooooooooooo 
♦ ♦*♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

III ill III 111 m HI LU iu II. HI m Hi UI III III m m II > III HI III ||| u i in 

O9>(n4®W0O0ini»lH0'lA(3HinNNHHp 
g)mcoO'cocp»oo'in- 9 r>«^mr>H 4 ^>pgmprtr*>g 
OtPip(>iAlA41A4N4«MANinMMno 40c 
0IAPl(P(P(P(PlP(PNNHNIA4P'lfl0440ll 
OOOOOOOOOOOOOO'O'O'OO'IA-HCOO 


<4 94 
AJ O' 
O' o 
43 O' 
r- 43 


oooooooo 
oooooooo 
+♦♦♦♦♦♦♦ 
ii ■ uj Uj IU UJ UJ UJ IU 
IPNOMAIPIA 4 
N(P0404n4 
4 « 4 WNNIAlP 
OHHNPimwtP 
0»4Nh OOOO 


O' pg 
ia m 
co m 
m cn 
o o 


>00000000 

>00000000 

I UJ UJ UJ UJ UJ UJ uj UJ 
fflOIMJHJOH 
f>i/\43-nr-_r->CT 

p»o*0'0'4min4 

(P(ni*lNlPNiP4 

OOOOOOOO* 


o o 

NNNNNNNNNNNHHHHHHNNNHfPHHNNNNNNNNNNNNNNN 


p fp Ip fp fp (p 0 

0000000 

I I I I I I I 

U UJ UJ UJ IU LU IU 
ON 4 A H 44 
A 4 N (PU3 Nf* 
h r* 4 — N 4 « 

43 n* o o> 43 in r- 
h O P* ffl -44 


I I 


• I I ( » ♦ » 


rg pg rg 

000 

1 I I 


nj 0 pg < 
o o o < 

lit 


INN 
>00 
I I 


uj uj uj i 

miAN»IA4|pn4HlA4 4NIANH'4(PACONm4lPIA40<l 
a>r~ 43 m r- r- n ia ia p* >h ia ia on N4 — piaipiah44®p'NC 
*yr^or~'4pgeo*4momr--9u>cn-HpgocDco®»Hknrgp^r*— <o*o 
OffliPiA<ON444<OOiPh44N««H«4IAiAHffliPiP<Pf | 

co-H 4 >or'- 9 <T 3 r>o- 9 or-. 4 p'* 4 »n> 40 « 4 m<t)' 0 >j-AjrgrgpgA<r 


I I I I I I I I I I I I I I 


-r-or~43f»o<T' 
) m |A 4 P ffl N4 
*4 AHOffl 44 
j 4 N 4 H N N 4 
JO-fflHHPH 


<n ro 

0 o 

1 I 

UJ UJ 
-« p- 

o o 
*4 co 
<A >4 
O 43 


pg IM * 

0 O < 

1 I 
UJ UJ u 
041 
O CO r 

r- rg r 
S 4 > 


> m o 

J 43 4 > 

» *4 4 


.OO 

r~oorg*g,4r\j4-m— •mm<t)'4mr--cinO'-H—i-4>H-g-H— •-9-H-H«*0'-*rg-4in<nfn4'rg~Hg3 


9 rg rg <n p 
> o o o c 
I I I I 

J UI IU UJ u 
4> 43 P*- 

o* <0 -o 
*4 *n m 
o» «» 0* 
o* o* u> 


0 41 
h o> r 
^ *4 • 
ia 4 r 
43 to r 


>n nn 
>000 
ill 
J uj UJ UJ 
' rg rn 43 
■rtp*h 
3 pg m oj 
3 «o pg .4 
■OOO 


NNP4N — 440444 — 


Nm p p p m 

000000 

llllll 

IU UJ UJ UJ UJ UJ 

r* 4 ® ffl r* n 
<0 O' — <0 O -9 
«4 O* — • — * 4 *0 
^ rg O' rg rg co 
O P n id 4 4 
• • • • • • 
-44® 4- 


rg<n- 4 in<or~<DO>OH*rgm.*in«r*a>< 40 - 4 Ajin' 4 iAr»<DO>o-H 
hhhhhhhhNNNNnNnNNNPPPPPP *4 *4 


Nm4iA44'ffl40-Arm<4lA 

HHH--J4— — NAINAjrgN 


IAIAIAIAIAIAIAIAIAIAIAIAIAIAiaiaiAIAIAIAI 


pgi 93 * 4 iA«^cDO'OH 4 Ajm^us«r-oO' 0 - 4 rgm^«A 4 >r‘ 0 ' 0 - 4 pgm«m'Of M 09 ' 0 - 9 rgm^iA<AP M coo>o-*njin'Ow>or^o 

HHHHHHHHNNNNNNNNNNPPPnnnAn -*9-99- — NNNNNNAINNNPPrtWPrtWiAH 


NlA4IA'0f*B4O-NP4IA 

— — — — — — — — nnnn n rg 


r'r^p-p-r^r^r^r-P“P*r'f»P fc f-r»r“p»r- 4 “p»P»p*^-r*P“r-®cD®cDa 3 oDa 3 <ocDococfflcoa 3 a 3 coa>eo®oeoeococDeo«a>®cooo* 0 ' 0 * 0 ‘« 40 > 0 ‘ 0 *o* 0 * 0 ' 0 ' 0 ' 0 * 0 ' 


82 


Figure 13.- Continued 
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Figure 13.- Continued 
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Figure 13.- Continued 
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Figure 13.- Continued 
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Figure 13.- Concluded 








Measurement 


Vortex Velocity, 
3 components 


Vortex Temperature 


Distance of Probe 
behind T-33 

Vertical Displace- 
ment of Wake 

T-33 Weight . 

T-33 Lift 
Constancy 

T-33 Air Speed 


Air Density 


Air Turbulence 
Levels 


Air Humidity 

Air Temperature 

Vertical Air 
Currents 


Crow Instability 


Lear jet - Wake 
Interaction 

Vortex Dissipation 
Mechanisms 

tear jet Trajectory 


Instrument 


Comments 


Page 
Ref s . 


DISA 3-wire anemometer Should be calibrated both on 

ground and in flight 
Tape recorder speed should re- 
flect digitizing required 
Recommend using an oscilloscope 
with a "memory" 


19-21 

21-22 

21 


DISA Single-Wire 
anemometer 


Must have constant current elec- 
tronics for frequency response 20 


LRC airborne DME 


22 


Altimeters 

Calculated . 

Pilot assessment 

Cockpit instruments 


Difference the Learjet and T-33 
altimeter readings 22 

Required to get theoretical T o 18 

Vertical acceleration of T-33 
should be minimized 18 

Value and constancy should be 
monitored and tape recorded 18 


Learjet ambient temp Required to get theoretical r Q 
and pressure probes • 18 


DISA 3-wire 
anemometer 
DISA Voltmeter 
Vanes on nose boom 


Problem is low output voltage 

13-15 

Has a 1 millivolt scale 14-15 

Limited to about 15 Hz 15-16 


On-Board meter or Ambient levels and within vortex, 

sampler if possible 25,31-36 


Learjet Temperature Need temperature distribution 

probe over ~ 2000 feet height 24 


Learjet Accelerometers Any estimates of vertical air cur- 
and hot-wire probes rents will be done from indirect 

measurements or isotropy of 
turbulence 24-25 


Movie camera 


Fly Learjet along wake and take 
movies 23 


Movie Camera 

Movie camera and 
Learjet penetration 


Done from chase plane filming the 
Learjet penetration 24 

Film dissipation region and fly 
Learjet probes through region 24 


Accelerometers on Single and double integration 

board Learjet necessary to establish position 

versus time for wake probes and 
corrections to hot-wire velocities — 


Figure 16.- Summary of recommendations on data collection. 



